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Abstract 


Shared-memory  multiprocessors  commonly  use  shared  variables  for  synchronization.  Our 
simulations  of  real  parallel  applications  show  that  large-scale  cache-coherent 
multiprocessors  suffer  significant  amounts  of  invalidation  traffic  due  to  synchronization. 
Large  multiprocessors  that  do  not  cache  synchronization  variables  are  often  more  severely 
impacted.  If  this  synchronization  traffic  is  not  reduced  or  managed  adequately, 
synchronization  references  can  cause  sever  congestion  in  the  network.  We  propose  a  class 
of  adaptive  backoff  methods  that  do  not  use  any  extra  hardware  and  can  significantly 
reduce  the  memoiy  traffic  to  synchronization  variables.  These  methods  use 
synchronization  state  to  reduce  polling  of  synchronization  variables.  Our  simulations  show 
that  when  the  number  of  processors  participating  in  a  barrier  synchronization  is  small 
compared  to  the  time  of  arrival  of  the  processors,  reductions  of  20  percent  to  over  95 
percent  in  synchronization  traffic  can  be  achieved  at  no  extra  cost.  In  other  situations 
adaptive  backoff  techniques  result  in  a  tradeoff  between  reduced  network  accesses  and 
increased  processor  idle  time. 
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Abstract  rr\ 

A  two-dimensional  array  is  called  monotone  if  the  maximum  entry  in  its  i-th  row 
lies  directly  below  or  to  the  right  of  the  maximum  entry  in  its  (i  -  l)-st  row.  (If  a  row 
has  several  maxima,  then  we  take  the  leftmost  one.)  A  two-dimensional  array  is  called 
totally  monotone  if  every  2x2  subarray  (t.e.,  every  2x2  minor)  is  monotone.  Totally 
monotone  arrays  were  introduced  by  Aggarwal.  Klawe,  Moran,  Shor,  and  Wilber,  who 
showed  that  several  problems  in  computational  geometry  could  be  reduced  to  the  prob¬ 
lem  of  finding  the  maximum  entry  in  each  row  of  a  totally  monotone  array.  They  also 
gave  a  sequential  algorithm  for  computing  the  row  maxima  of  an  n  x  n  totally  monotone 
array  in  0(n)  time.  In  this  paper,  we  generalize  the  notion  of  two-dimensional  totally 
monotone  arrays  to  multidimensional  arrays,  present  sequential  algorithms  for  finding 
maxima  in  such  arrays,  and  exhibit  a  wide  variety  of  problems  (involving  computa¬ 
tional  geometry,  dynamic  programming,  and  VLSI  river  routing)  that  can  be  solved 
efficiently  using  these  array-searching  algorithms. 


1  Introduction 


1.1  Motivation 

A  two-dimensional  array  A  =  {a*,,}  is  called  monotone  if  the  maximum  entry  in  its  i-th  row 
lies  directly  below  or  to  the  right  of  the  maximum  entry  in  its  (i  —  l)-st  row.  (If  a  row  has 
several  maxima,  then  we  take  the  leftmost  one.)  A  is  called  totally  monotone  if  every  2x2 
subarray  (i.e.,  every  2x2  minor)  is  monotone.  In  other  words,  for  all  t  <  k  and  j  <  l, 
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Figure  1.1:  For  i  <  k  and  j  <  l ,  d (pi,q:)  +  d (pk,qi)  >  d(p,,qi)  +  d (p*,g>). 

a,  j  <  a,,/  implies  <  ajc.i-  The  roin  maxima  problem  for  A  is  that  of  finding  the  maximum 
entry  in  each  row  of  the  array.  Now  suppose  some  oracle  has  an  n  x  m  totally  monotone 
array  A  and  that  we  wish  to  solve  the  row  maxima  problem  for  A  by  asking  the  oracle  for 
as  few  entries  of  A  as  possible.  Aggarwal  et  al.  [AKM*87]  showed  that  0(n  +  m)  queries  to 
the  oracle  suffice. 

Although  the  problem  of  finding  the  row  maxima  in  a  two-dimensional  totally  monotone 
array  may  seem  rather  odd  at  first  glance,  [AKM*87]  showed  that  several  problems  in  com¬ 
putational  geometry  can  be  reduced  to  one  or  more  instances  of  this  problem.  The  following 
example  illustrates  one  such  reduction  (borrowed  from  [AKM*87]). 

Suppose  we  are  given  a  convex  polygon  and  that  we  divide  it  into  two  convex  chains 
P  and  Q  (containing  n  and  m  vertices,  respectively)  by  removing  two  edges,  as  is  shown 
in  Figure  1.1.  Let  pi,...,pn  denote  the  vertices  of  P  in  counterclockwise  order  and  let 
qu .. .  ,qm  denote  the  vertices  of  Q  in  counterclockwise  order.  Then  for  1  <  i  <  k  <  n  and 
1  <  j '<  l  <  m,  consider  the  quadrilateral  formed  by  pi,  pk ,  9j,  and  qi.  By  the  quadrangle 
inequality  (which  states  that  the  sum  of  the  lengths  of  the  diagonals  of  any  quadrilateral  is 
greater  than  the  sum  of  the  lengths  of  any  pair  of  opposite  sides),  we  have 

d(P«,<Z;)  +  d(Pk>9!)  >  d(pM?<)+  '(:<  9j)- 

Thus,  if  we  imagine  an  n  x  m  array  A  =  {atiJ}  where  o,-j  -he  Euclidean  distance  from 
vertex  p,  €  P  to  vertex  q}  €  Q,  then  this  array  is  totally  monotone,  since  we  cannot  have 
both  ajj  <  eq,i  and  ak,j  >  ak,i.  Also,  since  any  entry  of  this  array  can  be  computed  in 
constant  time  (it  being  the  Euclidean  distance  between  two  points),  asking  for  an  entry  from 
an  oracle  simply  corresponds  to  evaluating  the  distance  between  the  corresponding  points. 
Hence,  by  using  [AKM*87],  we  can  find  the  farthest  vertex  in  Q  for  every  vertex  in  P  by 
evaluating  only  0(n  +  m)  distances.  In  fact,  [AKM*87]  showed  that  the  time  required  in 
addition  to  that  for  evaluating  these  0(n  +  m)  distances  is  also  linear  in  n-fm.  This  implies 
that  the  farthest  neighbor  problem  for  convex  chains  can  be  solved  in  linear  time  and  solves 
an  open  problem  in  [Tou83j. 
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1.2  Previous  Results 

As  the  previous  example  illustrates,  the  notion  of  total  monotonicity  is  closely  related  to 
the  quadrangle  inequality.  (It  is  this  relation  that  makes  the  combinatorial  formulation  of 
problems  in  terms  of  totally  monotone  arrays  such  a  useful  one.)  Another  related  property 
of  two-dimensional  arrays,  even  more  closely  tied  to  the  quadrangle  inequality,  was  studied 
previously  in  the  context  of  certain  transportation  problems.  In  the  late  eighteenth  century, 
G.  Monge  [MonSl]  observed  that  if  unit  quantities  (cannonballs,  for  example)  need  to  be 
transported  from  locations  X  and  Y  (supply  depots)  in  the  plane  to  locations  Z  and  W 
(artillery  batteries),  not  necessarily  respectively,  in  such  a  way  as  to  minimize  the  total 
distance  traveled,  then  the  paths  followed  in  transporting  these  quantities  must  not  intersect. 
In  1961,  A.J.  Hoffman  [Hof61]  elaborated  upon  this  idea  by  calling  an  n  x  m  array  C  =  {c^} 
a  Monge  array  if  c,i}  +  cl+]j+i  <  c,J+j  -f  Ci+\,j  for  1  <  i  <  n  and  1  <  j  <  m  and  by  providing 
a  greedy  algorithm  for  the  transportation  problem  when  the  cost  array  C  =  {ci,j}  is  a 
Monge  array1.  Monge  and  Hoffman’s  observations  together  imply  a  greedy  algorithm  for  the 
transportation  problem  when  the  sources  lie  on  one  line  in  the  plane  and  the  sinks  lie  on  a 
second,  parallel  line. 

Note  that  an  equivalent  way  to  define  an  n  x  m  Monge  array  C  —  {c,j}  is  to  require 
Ci,j  +  Ck,i  <  c,,;  +  Cfc,j  for  1  <  i  <  k  <  n  and  1  <  j  <  l  <  m,  as  the  following  proposition 
shows. 

Proposition  1.1  If  C  —  {c,,j}  is  an  n  x  m  Monge  array ,  then  for  1  <  i  <  k  <  n  and 
1  <  j  <  l  <  m,  Cij  +  cki,  <  a, i  +  ckj. 

Proof  Consider  any  i,  j,  k,  and  l  such  that  1  <  i  <  k  <  n  and  1  <  j  <  l  <  m.  By  the 
definition  of  a  Monge  array, 


ca,t  +  c*+i,t+i  <  c,.<+1  +  c4+iif 
for  1  <  s  <  n  and  1  <  t  <  m.  Thus,  for  1  <  t  <  m. 


fc-i  *-i 

+  cJ+i,t+i)  <  YAc.,+l  -I-  C4+i,t). 

f=»  1=1 

Cancelling  identical  terms  from  both  sides  of  this  inequality,  we  obtain 

c»',(  +  ck,t+ 1  ^  ^.t+l  +  C/i'f 

Consequently, 

i-i  i-i 

+ ct.t+i)  <  y^(cj,t+i + cfc.t). 

t=j  <=j 

Again  cancelling  identical  terms,  we  obtain 


Ci,j  "f"  Cfc.i  £  C»,l  "1"  Ck,y> 

‘In  fact,  Hoffman  [HofBl]  provided  a  necessary  and  sufficient  condition  that  the  cost  array  C  must  satisfy 
for  his  greedy  algorithm  to  work.  Furthermore,  Alon,  Cosares,  Hochbaum,  and  Shamir  [ACHS88]  have 
recently  developed  an  efficient  algorithm  for  testing  whether  this  necessary  and  sufficient  condition  holds. 
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the  desired  result.  ■ 


We  will  follow  Hoffman’s  terminology  and  call  c, j  -f  c*t/  <  c,./  +  the  Monge  condition 
and  Cij  ■+■  Ckj  >  Cij  +  Ckj  the  inverse  Monge  condition.  Note  that  the  inverse  Monge  condition 
implies  total  monotonicity,  but  not  vice-versa.  Although  total  monotonicity  is  a  somewhat 
weaker  constraint  than  the  inverse  Monge  condition,  it  turns  out  that  all  of  the  applications 
given  in  this  paper,  as  well  as  those  given  in  [AKM*87,  Wil88],  obey  the  Monge  or  inverse 
Monge  condition.  However,  we  usually  require  only  total  monotonicity  to  obtain  our  results. 

We  will  call  an  array  Monge  if  it  follows  either  the  Monge  condition  or  the  inverse  Monge 
condition,  since  an  array  satisfying  the  Monge  condition  may  be  converted  into  an  array 
satisfying  the  inverse  Monge  condition  (and  vice-versa)  by  reversing  the  order  of  its  columns. 
Note  that  in  a  similar  fashion,  we  can  flip  the  signs  in  our  definition  of  total  monotonicity 
(at  least  in  the  case  of  two-dimensional  arrays)  and  convert  back  and  forth  between  totally 
monotone  arrays  and  arrays  A  =  {a;,;}  such  that  a,j  >  a,,i  implies  ak,j  >  a.k,i  for  all  i  <  k 
and  j  <  l.  Furthermore,  the  problem  of  finding  the  maximum  entry  in  each  row  of  a  two- 
dimensional  array  A  and  the  problem  of  finding  the  minimum  entry  in  each  row  of  A  (the 
row  minima  problem)  are  equivalent  —  we  can  convert  one  to  the  other  by  simply  negating 
the  entries  of  A  and  reversing  the  order  of  its  columns.  We  make  use  of  these  dualities 
throughout  this  paper. 

As  mentioned  earlier,  totally  monotone  arrays  were  first  introduced  by  Aggarwal,  Klawe, 
Moran,  Shor,  and  Wilber  in  [AKM*87],  who  provided  a  linear  time  algorithm  for  com¬ 
puting  the  row  maxima  in  such  arrays.  Previous  applications  of  this  array- searching  algo¬ 
rithm  include  the  following.  [AKM*87]  showed  that  some  channel  routing  problems  can  be 
solved  in  linear  time  using  these  techniques,  thereby  improving  previous  0(n  lgn)  time  algo¬ 
rithms.  Aggarwal  and  Suri  [AS87]  used  array-searching  to  find  the  largest  empty  rectangle 
in  0(nlg2n)  time.  Wilber  [Wil88]  applied  array- searching  techniques  to  an  instance  of  dy¬ 
namic  programming  and  showed  that  the  concave  least-weight  subsequence  problem  can  be 
solved  in  linear  time;  this  improved  the  result  of  Hirschberg  and  Larmore  [HL87].  Finally, 
Klawe  and  Kleitman  [KK88]  used  array-searching  to  improve  a  previous  result  of  [AK88]  in 
computational  geometry. 

1.3  Main  Results  of  this  Paper 

In  this  paper,  we  present  a  framework  that  allows  us  to  efficiently  solve  a  wide  variety  of 
problems  involving  the  quadrangle  inequality.  We  limit  ourselves  to  sequential  computation; 
results  in  the  realm  of  parallel  computation  are  given  in  [AP89J.  This  paper  derives  its 
primary  inspiration  from  [AKM*87]  and  [YaoSO];  it  generalizes  and  incorporates  several 
results  provided  in  these  papers.  (The  relation  of  [AKM*87]  to  this  paper  will  be  observed 
throughout,  but  that  of  [YaoSO]  is  more  implicit.  We  assume  total  familiarity  of  the  reader 
with  [AKM*87j  and  (Yao80j.) 

The  balance  of  this  paper  is  organized  as  follows. 

Section  2  introduces  the  notion  of  multidimensional  monotone  arrays  and  provides  se¬ 
quential  algorithms  for  searching  in  such  arrays.  We  use  the  problem  of  finding  a  maximum 
perimeter  d- gon  inscribed  in  a  given  convex  n-gon  as  a  prototype  example.  For  this  problem, 
we  achieve  the  time  bound  of  [AKM*87]. 
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Section  3  applies  the  monotone  array  framework  to  the  problem  of  finding  minimum  area 
and  minimum  perimeter  circumscribing  polygons.  For  the  minimum  area  circumscribing 
d- gon  problem,  we  obtain  an  0(dn  4-  nlgn)  time  algorithm,  thereby  improving  the  best 
previous  result  of  [AKM*87]  by  a  factor  of  0((n  lg  d)/(d  +  lgn)).  For  the  minimum  perimeter 
circumscribing  triangle  problem,  we  obtain  an  0(n  lg  n)  time  algorithm,  thereby  improving 
the  result  of  [DePS7]  by  an  0(n2/  lgn)  factor. 

Section  4  presents  efficient  algorithms  for  dynamic  programming  problems  that  obey  the 
Monge  condition  or  the  inverse  Monge  condition.  We  give  an  0(n 2)  time  algorithm  for 
Frances  Yao’s  [YaoSO]  dynamic  programming  problem  using  the  monotone  array  framework. 
Although  our  algorithm  is  no  better  than  Yao’s  in  terms  of  asymptotic  complexity  (and 
probably  worse  in  practice),  it  provides  insight  into  other  dynamic  programming  problems 
that  obey  Monge  conditions.  We  then  use  this  insight  to  reduce  the  time  complexity  of  a 
particular  dynamic  programming  problem  related  to  biology  from  0(n2  lg2  n)  [EGG88]  to 
0(n2  lgn). 

Section  5  presents  efficient  sequential  algorithms  for  river  routing  in  VLSI.  We  show  that 
the  results  of  [AKM*87]  can  be  used  to  solve  the  offset  range  problem  in  linear  time  under 
very  weak  assumptions  regarding  the  routing  rules  for  wires.  This  generalizes  a  recent  result 
of  Siegel  [SieSS],  who  was  able  to  obtain  a  linear  time  algorithm  for  only  a  very  restricted 
class  of  wiring  rules. 

Section  6  discusses  open  problems. 


2  Multidimensional  Monotone  Arrays 

2.1  Basic  Definitions  and  Algorithms 

We  begin  by  extending  our  notions  of  monotone,  totally  monotone,  and  Monge  arrays  to 

higher  dimensions.  For  d  >  3,  let  A  =  {aiui2 . be  an  nj  x  n2  x  •  •  •  x  array.  Let 

*2(21  )•»  •  •  •  denote  the  second  through  d-th  coordinates  of  the  maximum  entry  in  the 

(d  —  1)- dimensional  plane  consisting  of  those  entries  whose  first  coordinate  is  *i,  i.e., 

a*i,»i(«j)--.*<i(«i)  =  -ma?  a»i 

(If  the  plane  corresponding  to  i2  contains  multiple  maxima,  we  chose  the  first  of  the  maxima 
ordered  lexicographically  by  their  second  through  d-th  coordinates.)  We  call  these  entries 
plane  maxima. 

Definition  2.1  A  is  monotone  if 

1.  for  1  <  ti  <  ji  <  ni,  we  have  t'*(t'i)  <  ik(j\)  for  all  k  between  2  and  d,  and 

2.  for  1  <  ij  <  «j,  the  (d  —  l) -dimensional  plane  of  A  consisting  of  those  entries  whose 
first  coordinate  is  ix  is  monotone. 

Definition  2.2  A  is  totally  monotone  if  every  2  x  2  x  ■  •  ■  x  2  d-dimensional  subarray  of  A 
is  monotone. 
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Definition  2.3  A  satisfies  the  Monge  condition  if  every  two-dimensional  plane  of  A  corre¬ 
sponding  to  fixed  values  of  d—  2  of  A's  d  coordinates  satisfies  the  Monge  condition.  Similarly, 
A  satisfies  the  inverse  Monge  condition  if  every  two-dimensional  plane  of  A  corresponding 
to  fixed  values  of  d  —  2  of  A's  d  coordinates  satisfies  the  inverse  Monge  condition.  Finally, 
A  is  Monge  if  it  satisfies  either  the  Monge  condition  or  the  inverse  Monge  condition. 

Note  that  every  subarray  of  a  totally  monotone  array  is  also  totally  monotone.  Similarly, 
every  subarray  of  an  array  satisfying  the  Monge  condition  satisfies  the  Monge  condition,  and 
every  subarray  of  an  array  satisfying  the  inverse  Monge  condition  satisfies  the  inverse  Monge 
condition. 

Proposition  2.4  Suppose  A  satisfies  the  inverse  Monge  condition,  and  consider  any  two 

entries  <z,lil2 . id  and  . J2 . u  of  A.  For  1  <  k  <  d,  let  xk  =  min{ik,jk}  and  let  yk  = 

max{zfc,jfc}.  Then 


axi  ,T2,...,Xd  d"  ay,  ,V2 yd  ^  ai\ ,ij, Oj,  ,j2l... ,jd  • 

Proof  We  prove  this  proposition  by  induction  on  d.  For  the  base  case  of  d  =  2,  the 
proposition  follows  immediately  from  the  definition  of  a  two-dimensional  array  satisfying  the 
inverse  Monge  condition. 

Now  suppose  the  proposition  holds  for  all  (d  —  l)-dimensional  arrays  satisfying  the  in¬ 
verse  Monge  condition,  and  consider  a  d-dimensional  array  A  satisfying  the  inverse  Monge 

condition  and  any  two  entries  a,,,i2t...,,d  and  a;iij2 . jd  from  A.  Without  loss  of  generality,  we 

assume  <  jx.  The  proof  then  breaks  down  into  two  cases. 

Case  1  For  all  k  between  2  and  d,  ik  >  jk  (i.e.,  xk  =  jk  and  yk  =  ik). 

Consider  the  (d  —  1)- dimensional  subarray  of  A  containing  those  entries  whose  second 
coordinate  is  i 2  and  the  (d  —  l)-dimensional  subarray  of  A  containing  those  entries  whose 
second  coordinate  is  j2.  Since  every  subarray  of  an  array  satisfying  the  inverse  Monge 
condition  satisfies  the  inverse  Monge  condition,  we  can  invoke  the  inductive  hypothesis  on 
these  subarrays  and  obtain 


‘»1.I2.J3  .  — .Jd 


+  a 


ji 


> 


+  a 


Jl  M<J  3, —,3d 


and 


+  a 


J1  ,tt.«3." 


,,  > 


a. 


»1  j2.«3.—.*if 


+  a 


ji  j3,}3,--o<r 


Similarly,  we  can  invoke  the  inductive  hypothesis  on  the  two-dimensional  subarray  containing 
those  entries  whose  third  through  d-th  coordinates  are  i3, . . . ,  ij,  respectively,  and  on  the 
two-dimensional  subarray  containing  those  entries  whose  third  through  d-th  coordinates  are 
j3, . . .  ,jj,  respectively.  This  gives  us 


..72  .'3. 


.Id 


+  a 


J)  .‘2  .»3. 


,d  > 


‘Il.l2.l3 . «d 


+  a 


}\  02.«3 . *d 


and 


**1  J2  03.«*^d 


+  0. 


Jl 


■Od  > 


‘il.i2.J3.— .Jd 


+  a 


Jl  .J 2  .J3 . Jd  ’ 
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Summing  these  four  inequalities  and  cancelling,  we  find 


“fl|l  .J2.J3 . 3d  "h  . id  —  “a'J.>2.I3 . id  "h  “GJl  ,33.33— ,3d' 

Since  Xj  =  tj,  j/i  =  j2.  and  for  2  <  k  <  d,  xk  =  jk  and  yk  =  t*,  this  gives  the  desired  result. 
Case  2  There  exists  an  /,  2  <  l  <  n,  such  that  t'/  <  jf  (t.e.,  x/  =  tj  and  yi  =  ji). 


Consider  the  (d—  l)-dimensional  subarray  containing  those  entries  whose  first  coordinate 
is  and  the  ( d  —  1  )-dimensional  subarray  containing  those  entries  whose  first  coordinate  is 
jj .  By  applying  the  inductive  hypothesis  to  these  subarrays,  we  obtain 


a. 


•  1.XJ.I3, 


,Id 


+  Uit 


V2.V3. 


;Vd  — 


X*1  .12  .*3 . >d 


+  a 


•1.33 .33 . J d 


and 

ajl,X2,X3 . Id  d-  ajl,V2.»3 . Vd  —  ajl.>2.'3 . «d  *h  °J1 ,33,33,— ,3dm 

Similarly,  by  applying  the  inductive  hypothesis  to  the  (d—  l)-dimensional  subarray  containing 
those  entries  whose  /- th  coordinate  is  and  to  the  ( d  —  l)-dimensional  subarray  containing 
those  entries  w’hose  /-th  coordinate  is  ji,  we  obtain 


1 1  .I3,r3' 


—id  d"  aji 


,•3,  <3 . id  — 


.>2,13. 


-f  a 


31, *3  ,X3,—.Id 


and 


ail, 33,33 . 3d  d"  ajl,V3.y3 . yd  —  °«».W.V3 . Vd  d"  a 31, 33, 33 3d' 


‘U.I2.X3 . Id  d"  3,y: . ltd  ^  -aU.<2.<3 . <d  d"  ^aji  ,33 ,33 ,— ,3d  ‘ 


Summing  these  four  inequalities  and  cancelling,  we  find 
2a, 

Since  Xj  =  ij  and  yx  =  j\ ,  this  gives  the  desired  result.  ■ 

Lemma  2.5  If  A  satisfies  the  inverse  Monge  condition,  then  A  is  totally  monotone. 


Proof  We  again  use  induction  on  d.  For  the  base  case  of  d  =  2,  the  lemma  follows 
immediately,  since  every  two-dimensional  array  satisfying  the  inverse  Monge  condition  is 
totally  monotone. 

Now  suppose  the  lemma  holds  for  all  (c/  —  l)-dimensional  arrays,  and  consider  a  d- 
dimensional  array  A  satisfying  the  inverse  Monge  condition.  Since  every  subarray  of  an 
array  satisfying  the  inverse  Monge  condition  satisfies  the  inverse  Monge  condition,  every 
(d—  l)-dimensional  plane  of  A  consisting  of  those  entries  whose  first  coordinate  is  some  fixed 
value  satisfies  the  inverse  Monge  condition.  Thus,  by  the  induction  hypothesis,  each  of  these 
planes  is  totally  monotone.  This  means  A  satisfies  the  second  property  of  totally  monotone 
arrays. 

Now  suppose  A  does  not  satisfy  the  first  property  of  totally  monotone  arrays,  t.e.,  there 
exist  i j,  and  /,  1  <  ii  <  <  nj  and  2  <  l  <  d,  such  that  ii(ii)  >  tj(ji).  For  2  <  k  <  d, 

let  xk  =  min{tfc(ti),t'fc(ji)}  and  let  yk  =  max{i*(*i), t*Ui)}.  By  Proposition  2.4, 


.13 . Id 


i 


This  gives  us  a 


ntradiction,  as  the  definition  of  plane  maxima  requires 


. >  G«l.*2 . Id 


and 


. —  °Jl,V2 . ild- 


Note  that  the  converse  of  Lemma  2.5  is  false,  i.e.,  totally  monotone  arrays  do  not  neces¬ 
sarily  satisfy  the  inverse  Monge  condition.  Thus,  total  monotonicity  is  a  weaker  condition, 
just  as  it  was  in  the  two-dimensional  case. 

We  call  the  problem  of  computing  the  plane  maxima  of  an  array  A  the  plane  maxima 
problem  for  A.  The  following  theorem  gives  an  upper  bound  on  the  time  necessary  to  solve 
the  plane  maxima  problem  for  a  totally  monotone  array. 

Theorem  2.6  For  d  >2,  the  plane  maxima  problem  for  an  x  n2  x  •  •  •  x  nd  d-dimensional 
totally  monotone  array  A  can  be  solved  in  0((n d  +  Ud_i)n*=i  lgnfc)  h'me. 

Proof  The  proof  is  by  induction  on  d.  For  the  base  case  of  d  =  2,  we  simply  apply  the 
O(n)  time  algorithm  of  [AKM*87j.  For  larger  values  of  d ,  we  use  a  simple  divide-and- 
conquer  approach.  We  begin  by  computing  the  maximum  entry  in  the  plane  corresponding 
to  i\  =  [ni/2j.  By  induction,  this  can  be  done  in  0((nj  +  nrf-i)rifc=2^gn*)  —  we 
simply  compute  the  plane  maxima  in  this  plane  and  then  spend  an  additional  0(n)  time 
to  compute  the  maximum  of  these  maxima.  This  gives  us  for  2  <  k  <  d.  Now 

since  A  is  monotone,  we  know  that  for  j,  <  i'i,  we  must  have  ik(ji)  <  tfcOi)  for  2  < 
k  <  d,  and  similarly  for  j\  >  tj,  we  must  have  ik(ji)  >  **(*i)  for  2  <  k  <  d.  Thus, 
we  need  only  consider  two  smaller  arrays,  one  (ij  —  1)  x  i2(m)  x  •••  x  id(i i)  and  the  other 
(ni  —  *i)  x  (n2  —  t2(*i)  + 1)  x  •  •  •  x  (nj  — i‘d(j i)  + 1),  in  computing  the  remaining  plane  maxima. 
This  gives  a  recurrence  for  the  time  T'(n1,n2, . . . ,  nd)  necessary  to  solve  the  plane  maxima 
problem  for  A  whose  solution  is  0((nd  -)-  n<j_i)  rifc=i  lg ”*)•  ■ 

Note  that  the  only  place  we  take  advantage  of  A’s  totally  monotonicity  is  in  solving  the 
row  maxima  problem  for  certain  two-dimensional  subarrays  —  we  can  do  almost  as  well  (i.e., 
compute  the  plane  maxima  of  A  in  0{nd  n^I2  lg  n*)  time)  using  only  monotonicity.  Also 
note  that  finding  the  plane  minima  in  a  totally  monotone  array  A  is  no  harder  that  finding 
the  plane  maxima  —  we  just  negate  all  of  A’s  entries  as  we  did  in  the  two-dimensional  case 
and  reverse  the  order  of  some  of  its  coordinates.  However,  we  cannot  convert  back  and  forth 
between  arrays  satisfying  the  Monge  and  inverse  Monge  conditions,  as  we  could  in  the  case 
of  two  dimensional  arrays. 

2.2  An  Example 

To  motivate  the  definitions  of  the  previous  subsection,  let  us  consider  a  concrete  example. 
The  maximum  perimeter  inscribed  d-gon  problem  is  defined  as  follows:  given  an  n-vertex 
convex  polygon  P,  find  a  d-gon  Q  contained  in  P  with  maximum  perimeter.  It  is  easy  to 
show  that  Q's  vertices  must  be  vertices  of  P.  Thus,  if  pi, . . .  ,pn  denote  the  vertices  of  P 
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and  per(?i,  t2. . . . ,  id)  denotes  the  perimeter  of  the  d-gon  corresponding  to  p,, , p,2 , . . .  ,  p,d,  we 

want  to  find  the  maximizing  per(ix,  ?2 _ ,  *<*)• 

For  the  case  of  d  =  3  (the  maximum  perimeter  inscribed  triangle  problem),  we  consider 
the  three-dimensional  array  A  =  where 

a  =|  »f  i  <  j  <k 

'.i.fc  —  oc  otherwise. 


If  we  can  find  the  maximum  entry  in  A ,  we  can  solve  the  maximum  perimeter  inscribed 
triangle  problem  for  P.  Note  that  we  do  not  explicitly  compute  all  of  the  entries  in  A; 
rather,  every  time  our  array-searching  algorithm  needs  a  particular  entry  from  A,  we  just 
calculate  the  perimeter  of  the  corresponding  triangle. 

Now  it  is  easy  to  verify  that  the  array  A  defined  above  satisfies  the  inverse  Monge 
condition.  Furthermore,  obtaining  a  particular  entry  in  the  array  requires  only  constant 
time.  Thus,  we  can  apply  Theorem  2.6  and  obtain  the  plane  maxima  of  A  in  0(n  lgrc)  time. 
We  can  then  compute  the  maximum  of  these  maxima  in  0(n)  additional  time,  w’hich  gives 
a  maximum  perimeter  inscribed  triangle. 

This  result  equals  the  result  obtained  by  Boyce  et  al.  in  [BDDG85];  moreover,  it  represents 
a  simpler  solution  to  the  problem,  as  Boyce  et  al.  require  a  number  of  additional  geometric 
properties  of  inscribed  triangles  that  complicate  their  proof. 

In  a  similar  fashion,  the  maximum  perimeter  inscribed  d-gon  problem  can  be  reduced  to 
the  problem  of  finding  the  maximum  entry  in  the  d-dimensional  array  A  =  . ,rf}  where 


f  per(?i,t‘2,  •  •  -Jd)  if  »i  <  t2  <  <  id 

(  —  oo  otherwise. 


It  is  easy  to  verify  that  A  satisfies  the  inverse  Monge  condition.  Furthermore,  0{d)  time  suf¬ 
fices  to  compute  any  entry  in  A.  Thus,  we  can  apply  Theorem  2.6  and  obtain  an  0(dn  lg^-2  n) 
time  algorithm  for  this  problem. 

It  remains  open  whether  one  can  do  better  than  O(nlgn)  time  for  the  inscribed  trian¬ 
gle  problem,  but  the  general  inscribed  d-gon  problem  can  be  solved  more  efficiently  if  we 
take  advantage  of  some  additional  structure  of  the  corresponding  d-dimensional  array;  this 
additional  structure  is  discussed  in  the  next  two  subsections. 


2.3  Monge-Composite  Arrays 

An  important  subclass  of  multidimensional  totally  monotone  arrays  consists  of  what  we  call 
Monge-composite  arrays.  As  one  might  expect,  an  array  is  Monge-composite  if  it  is  composed 
of  two-dimensional  Monge  arrays.  (These  arrays  may  be  either  variety  of  Monge  array,  ».e., 
they  can  obey  either  the  Monge  condition  or  the  inverse  Monge  condition,  but  they  all  must 
be  of  the  same  variety.)  More  precisely,  a  d-dimensional  Monge-composite  array  is  the  sum 
of  d-dimensional  extensions  of  two-dimensional  Monge  arrays,  where  we  define  the  sum  of 
two  arrays  and  the  extension  of  an  array  as  follows. 

Let  A  —  {a,, id)  and  B  =  {6<, be  nj  x  •••  x  nd  d-dimensional  arrays.  The  sum  of 

A  and  B  (written  A  +  B)  is  the  ni  x  ■  ■  •  x  nj  d-dimensional  array  C  =  {c„ I(<}  where 

C<1 . «d  =  °«j . id  d"  . id' 
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for  all  , id. 

Now  let  E  =  {etl . ,d}  be  an  nx  x  •  •  •  x  nd  d-dimensional  array.  For  any  dimension  k 

between  1  and  d  +  1  and  any  size  h,  the 

Hi  x  •  •  •  x  nk- 1  x  h  x  rik  x  •  •  •  x  nd 

(d  -f  l)-dimensional  array  k  =  {/,, . ,d+1}  is  an  extension  of  E  if 

/»]  .'*-1  .»*,»*+  1  ,.,*d+l  e'J  >*t+l  .....‘d+l  ’ 

for  all  i'i, . . .  ,id+i-  (E  is  just  n  copies  of  F,  each  one  a  plane  of  F  corresponding  to  some 
fixed  value  of  F’s  fc-th  coordinate.)  Furthermore,  any  extension  of  an  extension  of  E  is  also 
an  extension  of  E. 

From  these  definitions,  it  is  clear  that  each  entry  of  a  d-dimensional  Monge-composite 
array  A  —  may  be  written 


k<l 

where  for  all  k  <  /,  the  nk  x  nt  array  }  is  a  Monge  array. 

Proposition  2.7  The  sum  of  two  Monge  arrays  is  also  Monge. 

Proposition  2.8  The  two-dimensional  extension  of  a  vector  or  a  scalar  is  Monge. 
Proposition  2.9  Every  Monge-composite  array  A  is  Monge. 


Proof  We  must  show  that  all  two-dimensional  planes  of  A.  corresponding  to  fixed  values 
of  d  —  2  of  /4’s  d  coordinates,  are  Monge.  To  see  why  this  is  true,  consider  any  such  plane. 
This  plane  is  the  sum  of  a  two-dimensional  Monge  array,  some  vectors,  and  some  scalars; 
thus,  the  plane  is  Monge.  ■ 


Two  special  cases  of  Monge-composite  arrays  are  path-  and  cycle-decomposable  arrays. 
An  array  A  =  {a,, . ,d}  is  path-decomposable  if  each  of  its  entries  satisfies 


. «d 


=  wn  ,.2  +  + - bw, 


(d-l.d) 


where  the  re’s  are  entries  from  two-dimensional  Monge  arrays  as  before.  A  is  cycle- 

decomposable  if  each  of  its  entries  satisfies 


(1,2)  . 

G,, . =  U>,j  ,2  + 


.  (d-\,d)  ,  (d,l) 
4-  W.  ■  4-  10;  ’■  ’ . 

'  '  »d,*l 


Theorem  2.10  The  plane  maxima  of  an  n\  x  •  •  ■  x  nd  d-dimensional  path-decomposable 
array  A  can  be  computed  in  0(J2t=i  nk)  time. 
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Proof  The  proof  is  by  induction  on  d.  The  base  case  of  d  =  2  follows  immediately 
from  [AKM*87].  For  d  >  2.  we  assume  by  induction  that  the  theorem  holds  for  all  lower¬ 
dimensional  path-decomposable  arrays.  Since  A  =  {at] . is  path-decomposable,  we  can 

write 


■M 


(1.2)  .  (2.3)  , 

=  w  4*  w  ■  4 


•  +  w{d-lJ) 


where  the  w's  are  entries  from  two-dimensional  Monge  arrays.  Now  consider  the  n2  x  •  •  •  x  nrf 
(d  —  l)-dimensional  array  B  =  {6,2 . id)  where 


■'«! . *d 


(2.3)  . 
—  in'  '  4. 
>2.«3  ~ 


+  W 


(d-U) 


By  induction,  we  can  find  the  plane  maxima  of  B  in  0(£*_2n*)  time.  Since  the  maximum 
entry  in  the  plane  of  A  corresponding  to  a  particular  value  ?!  of  the  first  coordinate  is  just 


min 

«2 


(1.2)  , 

tu  +  mm 


(tc(2-3) 


+  ■  •  •  +  u>. 


•d-i  .‘d 


}}• 


we  need  only  find  the  row  maxima  a.  ihe  sum  of  and  the  appropriate  two-dimensional 

extension  of  the  vector  of  B  s  plane  maxima.  Since  this  sum  is  a  Monge  array,  we  can  find 
its  row  maxima  in  an  additional  0(nj  4-  n2)  time,  which  yields  the  desired  result.  ■ 


Theorem  2.11  The  plane  maxima  of  an  x  ■  ■  ■  nj  d-dimensional  cycle-decomposable  array 
A  can  be  computed  in  0((52*=2  n*)  Ig^i)  hme. 

Proof  First  note  that  each  plane  A-n  of  A,  corresponding  to  a  fixed  value  of  A's  first 
coordinate,  is  path-decomposable,  since 

Ah  =  H'(W)  +  +  •  •  •  + 

where  W7*2,3*  is  the  sum  of  W'(2,3*  and  two-dimensional  extensions  of  the  vectors  V 2  =  {in!1,2*] 

1  *J.*2  J 

and  Vd  =  {u^-11}. 

1  *d,*l  3 

This  means  we  can  compute  a  particular  plane’s  maximum  entry  in  0(]Cl_2  nk)  time. 
Combining  this  with  the  divide-and-conquer  approach  of  Theorem  2.6,  we  obtain  the  desired 
result.  ■ 


2.4  More  Examples 

Returning  to  th*  maximum  perimeter  inscribed  d-gon  problem,  it  is  easy  to  verify  that 
the  d-dimensie.)  array  A  defined  in  Subsection  2.2  is  cycle-decomposable.  Thus,  we  can 
apply  Theorer  1  and  obtain  an  0(dn  lgn)  time  algorithm  for  the  problem.  Furthermore, 
it  is  possible  to  reduce  this  time  bound  to  0(dn  -f  nlgn)  using  one  additional  technique 
from  [BDDG85]  v  c...  allows  us  to  restrict  our  attention  to  a  relatively  small  subarray  of 
A.  (We  preset*  this  technique  in  Subsection  3.1  in  the  context  of  the  minimum  perimeter 
circumscribing  d-gon  problem.) 

As  Boyce  et  al.  [BDDG85]  note,  the  same  approach  can  also  be  used  to  solve  the  max¬ 
imum  area  inscribed  d-gon  problem  in  0(dn  4- nlgn)  time.  Although  the  corresponding 
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d-dimensional  array  does  not  satisfy  the  inverse  Monge  condition  (nor  is  it  totally  mono¬ 
tone),  it  can  be  shown  that  we  only  need  to  consider  certain  subarrays  of  this  array,  and 
these  subarrays  do  satisfy  the  inverse  Monge  condition. 

Another  application  in  which  multidimensional  Monge-composite  arrays  naturally  arise  is 
the  following  problem:  given  a  convex  n-gon  P  with  vertices  Pi,P2,  -  ■  ■  ,Pn  in  counterclockwise 
order,  find  for  each  vertex  p,  of  P  the  longest  d-link  convex  path  from  pi  to  p,.  Corresponding 
to  this  problem  is  a  (d  —  l)-dimensional  path- decomposable  array,  whose  plane  maxima  are 
these  maximal  paths.  Thus,  we  can  apply  Theorem  2.10  and  obtain  these  paths  in  O(dn) 
time.  (This  result  is  implicit  in  [AKM*87].) 

Note  that  not  all  applications  involve  path-  or  cycle-decomposable  arrays;  consider,  for 
example,  the  problem  of  finding  a  maximum  weight  d- clique  of  vertices  from  a  convex  n-gon, 
where  the  weight  of  a  clique  is  defined  to  be  the  sum  over  all  pairs  of  vertices  in  the  clique 
of  the  distance  between  the  two  vertices.  An  instance  of  this  problem  corresponds  to  a  d- 
dimensional  Monge-composite  array,  each  of  whose  entries  may  be  computed  in  0(<P)  time, 
and  thus  may  be  solved  in  0(d2n  \gd~2  n)  time.  Furthermore,  the  d- dimensional  array  we 
consider  in  Subsection  3.2  in  connection  with  the  minimum  perimeter  circumscribing  triangle 
problem  is  not  even  Monge-composite. 


3  Finding  Minimum  Circumscribing  Polygons 

In  this  section,  we  apply  algorithms  for  searching  in  totally  monotone  arrays  to  the  problem 
of  finding  minimum  area  and  minimum  perimeter  circumscribing  polygons.  Given  an  n- 
vertex  convex  polygon  P  and  an  integer  d  between  3  and  n,  we  want  to  find  a  minimum 
area  or  minimum  perimeter  d-gon  Q  containing  P.  Note  that  for  both  area  and  perimeter 
minimization,  Q  must  clearly  be  convex,  and  each  of  its  d  edges  must  contact  P.  Also  note 
that  if  we  can  find  area-  and  perimeter-optimal  d-gons  circumscribing  convex  n-gons,  then  we 
can  also  find  area-  and  perimeter-optimal  convex  d-gons  containing  arbitrary  sets  of  points 
in  the  plane,  since  any  convex  polygon  containing  a  set  of  points  must  contain  the  points’ 
convex  hull. 

Chang  and  Yap  [CY84]  showed  that  a  minimum  area  circumscribing  d-gon  can  be  found  in 
0(n3lgd)  time  using  dynamic  programming.  Aggarwal,  Chang,  and  Yap  [ACY85]  improved 
this  result  to  0(n2lgdlgn)  time,  and  it  is  further  improved  to  <?(n2lgd)  in  [AKM*87]. 
We  extend  (in  a  non-trivial  manner)  the  techniques  of  Boyce  et  al.  [BDDG85]  for  finding 
maximum  perimeter  inscribed  d-gons  to  obtain  an  0(dn  -f-  nlgn)  time  algorithm  for  the 
minimum  area  circumscribing  d-gon  problem.  We  present  this  algorithm  in  Subsection  3.1. 

As  for  the  perimeter  minimization  problem,  the  only  known  results  are  for  d  =  3.  DePano 
[DeP87]  showed  that  the  minimum  perimeter  circumscribing  triangle  cam  be  computed  in 
0(n3)  time.  In  Subsection  3.2,  we  improve  this  to  O(nlgn)  time  using  a  straightforward 
application  of  our  techniques  for  searching  in  three-dimensional  totally  monotone  arrays. 

In  describing  our  algorithms  for  finding  minimum  area  and  minimum  perimeter  circum¬ 
scribing  d-gons,  we  use  the  following  conventions.  For  any  convex  m- gon  iZ,  we  let  vf*, . . . ,  v% 
and  ef , . . . ,  e*  denote  P's  vertices  and  edges,  respectively,  in  counterclockwise  order,  where 
ef  connects  vf1  and  u^+])modm.  We  use  the  letter  P  for  the  convex  n-gon  to  be  circumscribed 
and  the  letter  Q  for  convex  d-gons  circumscribing  P. 


12 


We  also  need  the  following  definitions. 

1.  If  an  edge  ef  of  a  circumscribing  polygon  Q  touches  P  (which  it  must,  if  Q  has  minimal 
area  or  perimeter),  then  its  contact  point  is  the  part  of  P  that  it  touches.  A  contact 
point  is  always  either  an  edge  or  a  vertex  of  P. 

2.  Two  circumscribing  polygons  <5  and  Q'  interleave  if  the  contact  points  of  Q  and  Q' 
alternate.  In  other  words,  between  every  two  consecutive  contact  points  of  Q  (in 
the  counterclockwise  ordering  of  vertices  and  edges  of  P)  is  a  contact  point  of  Q' 
(perhaps  one  of  the  two  consecutive  contact  points  of  Q).  Similarly,  between  every  two 
consecutive  contact  points  of  Q'  is  a  contact  point  of  Q. 

3.  An  edge  ef  of  Q  is  flush  with  P  if  its  contact  point  is  an  edge.  <5  itself  is  flush  with  P 
if  all  its  edges  are  flush  with  P. 

4.  An  edge  ef  of  Q  is  balanced  if  its  midpoint  lies  on  P. 

5.  An  edge  ef  of  Q  determines  two  half-planes;  let  Hi  denote  the  half-plane  that  does  not 
contain  P.  ef  is  a  c-edge  if  the  lines  containing  its  neighbors  e£.j  and  converge 
(t.e.,  intersect)  in  Ht  or  are  parallel;  otherwise,  ef  is  a  d-edge.  Equivalently,  ef  is  a 
c-edge  if  the  sum  of  the  two  internal  angles  of  Q  corresponding  to  ef ’s  endpoints  is 
greater  than  or  equal  to  7r,  and  ef  is  a  d-edge  if  this  sum  is  less  than  n. 


3.1  Area  Minimization 

Our  algorithm  for  finding  a  minimum  area  circumscribing  d- gon  has  two  parts.  First,  we 
restrict  our  attention  to  flush  circumscribing  d- gons  and  use  the  techniques  of  [BDDG85] 
to  find  one  with  minimal  area.  Then,  we  use  this  minimum  area  flush  d-gon  (and  a  lemma 
due  to  DePano  [DeP87])  to  obtain  a  circumscribing  d-gon,  possibly  not  flush,  whose  area  is 
minimal  among  all  circumscribing  d-gons. 


3.1.1  Finding  the  Best  Flush  d-gon 

The  techniques  given  by  Boyce  et  al.  [BDDG85]  for  finding  a  maximum  perimeter  inscribed  d- 
gon  can  also  be  used  to  find  a  minimum  area  flush  circumscribing  d-gon  in  0(dn  lg  n+n  Ig2  n) 
time.  (This  was  pointed  out  in  the  concluding  section  of  [BDDG85].)  Furthermore,  the 
techniques  of  [AKM*87]  reduce  the  time  complexity  of  this  problem  to  0{dn  +  nlgn).  For 
the  sake  of  completeness,  we  will  describe  this  result,  recasting  it  in  terms  of  multidimensional 
totally  monotone  arrays. 

For  1  <  t  <  n,  let  r?cw  be  the  ray  containing  ef  with  vf  as  its  origin,  and  let  rfw 
be  the  ray  containing  ef  with  as  its  origin.  (The  superscript  of  rfcw  indicates  that 
it  is  a  counterclockwise  “extension”  of  ef ,  and  the  superscript  of  r indicates  that  it  is  a 
clockwise  “extension.”)  If  r,ccw  intersects  r'w,  let  Rij  be  the  region  outside  P  bounded  by 
rfcw,  rjw,  and  the  edges  ef.], . . .  ,ef_j  of  P  (the  shaded  region  in  Figure  3.1).  Now  consider 
the  two-dimensional  arrays  W  =  {u\tJ}  and  W  =  {tuf,},  where 


area (Ri,j)  if  t  <  j  and  rfcw  intersects  r^w,  and 
oo  otherwise, 
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Figure  3.1:  If  if"  intersects  rc*,  then  Ritj  is  the  shaded  region  between  r,ccw,  rf\  and  P. 


Figure  3.2:  Suppose  u>,,,  in,/,  wkj,  and  in*/  all  correspond  to  regions  outside  of  P,  i.e.,  i  <  k  <  j  <  l 
and  r,ecw  intersects  rf*.  The  region  12,/  bounded  by  rf"  and  rf*  contains  both  the  region  ff,y 
bounded  by  r,cc*  and  rf*  (shaded  with  horizontal  lines)  and  the  region  Rki  bounded  by  rf*  and 
rf"  (shaded  with  vertical  lines).  Moreover,  the  intersection  of  Rij  and  if*/  is  exactly  the  region  if*j 
bounded  by  rf*  and  rf*. 


and 


j  area  (Rij)  if  i  >  j  and  r,cc*  intersects  rf*,  and 
1  oo  otherwise. 


Lemma  3.1  Both  W  and  W'  satisfy  the  Monge  condition. 


Proof  We  only  prove  the  lemma  for  W;  the  proof  for  W'  is  similar.  For  1  <  t  <  <  n 

and  1  <  j  <  l  <  «,  we  must  show  that  ujtJ  -f  t ukl  <  wu  4-  wkj.  We  consider  two  cases.  If 
either  wu  =  oo  or  wk:  =  oo,  then  the  Monge  condition  follows  immediately.  If,  on  the  other 
hand,  w«  and  wkj  both  correspond  to  regions  outside  of  P  (Rn  and  Rk],  respectively),  then 
Wjj  and  wki  must  also  correspond  to  regions  outside  of  P  ( Rij  and  Rki,  respectively),  and  the 
four  regions  must  overlap  as  in  Figure  3.2.  Now  consider  the  region  Rij  U  Rki-  It  has  area 
Wij  +  t vki  -  wkj,  since  the  intersection  of  R:  and  Rkt  is  exactly  Rkj.  Moreover,  R:  U  Rki  is 
contained  in  Ri,  which  implies  uf  +  wki  —  wkj  <  wj  or  Wij  +  u>*j  <  wj  +  wkj.  ■ 

By  summing  d-dimensional  extensions  of  W  and  W1 ,  we  obtain  the  d-dimensional  array 
A  =  {ailiia . id},  where 

=  d"  Wij.t's  +  •  ’  '  +  Wid-\,U  d" 
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A  is  clearly  a  Monge-composite  cycle-decomposable  array.  Furthermore,  A  contains  an 
entry  corresponding  to  every  possible  flush  circumscribing  d- gon.  Specifically,  the  area  of 
the  flush  circumscribing  d- gon  with  contact  points  efi,  e£, . . . ,  e£,  ij  <  i2  <  •  ■  •  <  t‘d,  is 

Ojj  . . -|-area(P).  Moreover,  only  those  entries  corresponding  to  circumscribing  d-gons  are 

less  than  oc;  thus,  to  find  a  minimum  area  flush  circumscribing  d- gon,  we  need  only  find  a 
maximum  entry  in  A. 

By  applying  Proposition  2.11  directly,  we  can  find  this  entry  in  0(dn  lg  n)  time.  However, 
this  time  complexity  can  be  reduced  to  0(dn  +  nlgn)  using  a  theorem  concerning  the  area 
of  interleaving  d-gons.  (The  theorem  we  prove  is  actually  significantly  more  general  than 
we  need  it  to  be  in  obtaining  an  0(dn  -f  n  lgn)  time  algorithm  for  the  minimum  area  flush 
circumscribing  d- gon  problem,  but  we  will  use  it  again  later  in  this  section  in  the  context  of 
two  other  related  problems.) 

Before  we  can  prove  this  theorem,  however,  we  first  need  a  few  more  definitions.  Let  Qa 
and  Qk  denote  circumscribing  d-gons.  We  define  an  edge  exchange  operation  as  follows.  Let 
E  denote  the  union  of  Qa  and  Qb's  edges.  To  exchange  edges  between  Qa  and  Qb ,  we  select 
a  d-edge  subset  E'  of  E.  and  form  two  new  circumscribing  d-gons,  the  first  consisting  of  the 
edges  of  E'  (extended  or  shortened  as  necessary  to  form  a  circumscribing  d-gon)  and  the 
second  consisting  of  the  edges  of  E  —  E'. 

Now  let  Qa  and  Qb  denote  sets  of  circumscribing  d-gons.  We  will  say  that  Qa  and  Qb 
are  closed  under  edge  exchange  if  for  any  d-gon  Qa  €  Qa  and  any  d-gon  Qb  6  Qb ,  any  edge 
exchange  between  Qa  and  Qb  produces  a  d-gon  in  Qa  and  a  d-gon  in  Qb. 

Theorem  3.2  Suppose  two  sets  Qa  and  Qb  of  circumscribing  d-gons  are  closed  under  edge 
exchange.  Furthermore,  suppose  that  Qa  has  minimum  perimeter  among  d-gons  in  Qa  and 
that  Qb  has  minimum  perimeter  among  d-gons  in  Qb.  Then  Qa  and  Qb  interleave. 

Proof  Suppose  Qa  and  Qa  do  not  interleave.  Let  a\,  a2, . . . ,  a<f  be  the  contact  points  of 
Qa,  and  let  &i,  62, . . .  ,64  be  the  contact  points  of  Qb.  Since  Qa  and  Qb  do  not  interleave, 
there  exists  at  least  one  pair  (aj_j,Oi)  of  consecutive  contact  points  of  Qa  such  that  no 
contact  points  of  Qb  lie  between  a,_  1  and  c,  (inclusive)  in  the  counterclockwise  ordering  of 
P's  vertices  and  edges.  There  also  exists  at  least  one  pair  (bj,  fy+i)  of  consecutive  contact 
points  of  Qb  such  that  no  contact  points  of  Qa  lie  between  bj  and  bj+i  (inclusive).  Moreover, 
there  exists  such  a  pair  (a<_i,at)  and  such  a  pair  (bj,  6J+1)  separately  only  by  alternating 
contact  points.  In  other  words,  there  exist  i,  j ,  and  k,  such  that  the  contact  points  between 
cii_i  and  6J+i  are  (in  order) 

li  b) fe)  ®i+ ^k+  li  0»'+2>  •  •  •  1  bj— 1>  bj,  bj+  j. 

,  - - — v«  ■  " 

alternating  contact  points 

Now  suppose  we  exchange  edges  between  Qa  and  Qb  and  form  a  d-gon  Q'a  with  contact 
points 

®i  1  •  •  ■  1  ®i-ii  bit,  bi 1 ,  ■  •  •  ,bj,  Oj-fc+»+ 1 1  •  •  •  1  Qd 

and  a  d-gon  Q'b  with  contact  points 

b  1 ,  .  .  .  ,  j ,  O^,  G,+i,  .  .  .  ,  Qj— Je+ti  bj+\ ,  •  •  •  ,  b,i. 
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Figure  3.3:  The  sum  of  the  areas  of  Qa  and  Qb  (indicated  by  the  dotted  lines)  is  exactly  the  sum 
of  the  areas  of  Q'a  and  Q[  (indicated  by  the  shaded  lines)  plus  the  area  of  the  two  shaded  regions. 

For  any  two  contact  points  c  and  c',  let  Rc%c>  denote  the  region  outside  P  bounded  by  P  and 
lines  through  the  edges  of  Qa  or  Qb  touching  c  and  c'.  As  is  suggested  in  Figure  3.3, 

area(Q'a)  +  area(Qi)  =  area(Q0)  +  area(Qb) 

—  area(Pai_li0l)  -  area (P„,_*+(,0,_*+l+1) 

—  area (Rbk_ubk)  -  area(f?Mj+I) 

+  area  (£„,_„**)  +  area  (Pb„a,_*+>+, ) 

+  area(f?tk_li0|)  +  area(Paj_k+tifcj+1). 

Since  bk_x  precedes  a,_i  in  the  counterclockwise  ordering  of  P’s  vertices  and  edges, 

area(Ptk_lia.)  +  area(Pai_1,6k)  <  area(P6k_,,6j  +  area(P0, 

by  the  same  argument  we  used  in  the  proof  of  Lemma  3.2.  Similarly, 

area (P0,_k+„6J+1 )  +  area(Pi)j,a;_k+(+1 )  <  area(PQ;_k+i,aj_k+i+1 )  +  area(P6„0,+J ). 

Thus,  area(Q^)  +  area(Q|,)  <  area(Qo)  +  area(Qb).  Since  Qa  and  £>&  are  closed  under  edge 
exchange,  one  of  the  new  d-gons  is  in  Qa  and  the  other  is  in  Qb.  Without  loss  of  generality, 
we  assume  Q'a  €  Qa  and  Q'b  €  Qb-  Now  either  area(<2'0)  <  area(@0)  or  area(Q^)  <  area(<5t), 
both  of  which  are  contradictions.  ■ 

Corollary  3.3  Let  Qi  be  a  minimum  area  flush  circumscribing  d-gon  with  ef  as  a  contact 
point.  Every  minimum  area  flush  circumscribing  d-gon  interleaves  Qi. 

Proof  Any  edge  exchange  between  a  circumscribing  d-gon  with  ef  as  a,  contact  point  and 
an  arbitrary  circumscribing  d-gon  produces  a  circumscribing  d-gon  with  ef  as  a  contact  point 
and  an  arbitrary  circumscribing  d-gon.  Thus,  this  corollary  follows  from  Thereom  3.2.  ■ 

Returning  to  the  problem  of  finding  a  minimum  area  flush  circumscribing  d-gon,  note 
that  finding  a  minimum  area  flush  circumscribing  d-gon  Qi  with  ef  as  a  contact  point  is 


equivalent  to  finding  a  minimum  entry  in  first  plane  of  A.  This  can  be  done  in  0{dn)  time 
(since  this  plane  is  path-decomposable).  Let  e£,...,e£  be  the  contact  points  of  Q\  (by 
definition,  j\  =  1;  by  convention.  jd+\  =  n).  These  edges  define  d  intervals  Jj, . . . ,  7j  of  edges 
from  P,  where  Ik  =  [e£,  e£  ].  Corollary  3.3  tells  us  that  we  need  only  consider  edges  in  Ik 
for  the  fc-th  contact  point  of  a  minimum  area  circumscribing  d-gon.  In  other  words,  we  need 
only  search  the  nj  x  n2  x  ■  •  •  x  n<j  subarray  of  A ,  nk  =  jk+i  —jk  +  1,  containing  those  entries 
a^i, . .  where  jk  <  4  <  jk+i  for  all  k  between  1  and  d.  Since 

En‘  =  °(n) 

k=l 

and  every  subarray  of  a  cycle-decomposable  array  is  also  cycle-decomposable,  we  can  use 
Theorem  2.11  to  find  a  minimum  entry  in  this  subarray,  corresponding  to  a  minimum  area 
flush  circumscribing  d-gon,  in  0{n\gn)  additional  time.  This  gives  the  entire  algorithm  a 
time  complexity  of  0(dn  +  n  Ig  n). 

3.1.2  Using  the  Best  Flush  d-gon  to  Obtain  the  Best  Arbitrary  d-gon 

In  [DeP87],  DePano  provides  the  following  geometric  characterization  of  minimum  area  cir¬ 
cumscribing  d-gons. 

Lemma  3.4  ([DeP87])  Let  P  be  any  convex  n-gon.  For  3  <  d  <  n,  if  Q  is  a  minimum 
area  d-gon  Q  circumscribing  P,  then  either 

1.  all  d  edges  of  Q  are  flush  with  P,  or 

2.  d  —  1  edges  of  Q  are  flush  with  P,  and  the  non-flush  edge  is  a  balanced  d-edge. 

This  lemma  allows  us  to  relate  minimum  area  flush  circumscribing  d-gons  and  mini¬ 
mum  area  arbitrary  circumscribing  d-gons.  Specifically,  we  have  the  following  corollary  to 
Theorem  3.2. 

Corollary  3.5  Let  Q'  be  a  minimum  area  flush  circumscribing  d-gon.  Every  minimum  area 
circumscribing  d-gon  Q  interleaves  Q' . 

Proof  Let  Q!  denote  the  set  of  all  flush  circumscribing  d-gons,  and  let  Q  denote  the  set 
of  all  circumscribing  d-gons  whose  first  d  —  1  edges  are  flush  with  P  and  whose  d-th  edge  is 
a  balanced  d-edge.  By  Lemma  3.4,  every  minimum  area  circumscribing  d-gon  is  in  Q  U  Q! . 
Moreover,  every  minimum  area  circumscribing  d-gon  has  minimal  area  among  d-gons  in 

<2u<2'. 

Now,  Q!  and  QuQ'  are  clearly  closed  under  edge  exchange.  Thus,  by  Theorem  3.2,  every 
minimum  area  circumscribing  d-gon  must  interleave  every  minimum  area  flush  circumscribing 
d-gon.  ■ 

Now  suppose  we  have  found  a  minimum  area  flush  circumscribing  d-gon  Q' ,  using  the 
techniques  of  [BDDG85].  Let  e£, . . . ,  e£  be  the  contact  points  of  this  d-gon.  Without  loss  of 
generality,  assume  j\  =  1  (if  it  is  not,  we  can  renumber  the  edges  of  P).  Also,  for  notational 
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convenience,  let  jd+i  =  n.  The  edges  ef , . . . ,  ef  define  d  intervals  of  edges  from 

P,  where  lk  =  [ef  ,ef+J.  Lemma  3.5  tells  us  that  we  need  only  consider  edges  or  vertices 
in  /jt  for  the  k- th  contact  point  of  a  minimum  area  circumscribing  d- gon  Q.  (A  vertex  is 
considered  to  be  in  Ik  if  both  the  edges  incident  to  it  are  in  Ik)  Furthermore,  we  can  use 
Lemma  3.4  to  show: 

Proposition  3.6  There  are  at  most  three  intervals  that  might  contain  the  contact  point  of 
the  non-flush  edge  of  a  minimum  area  circumscribing  d-gon  Q.  Moreover,  we  can  identify 
these  intervals  in  0(d)  time. 

Proof  For  each  edge  ef  of  P,  let  r?°*  be  the  ray  containing  ef  with  vf  as  its  origin.  Now 
for  each  interval  Ik ,  consider  the  ray  rf*  associated  with  the  interval’s  first  edge  ef .  Let  ak 
be  rf^’s  angle  with  respect  to  rjcw,  measured  in  a  counterclockwise  direction.  Clearly, 

0  =  Qj  <  02  <  •  •  •  <  ctd  <  2r. 

Now  suppose  the  non-flush  edge  of  Q  is  e^.  It  must  contain  an  edge  or  vertex  of  P  in 
the  interval  Ik.  Moreover,  the  contact  point  of  is  an  edge  ef  €  h-\  and  the  contact 
point  of  e*+1  is  an  edge  ef  €  h+ 1-  Since  ef  lies  in  /*_i,  the  angle  forms  with  rf0"  is  at 
least  o*_ i.  Similarly,  since  ef  lies  in  lk+ 1,  the  angle  r*,cw  forms  with  rf™  is  at  most  a*+2. 
Furthermore,  DePano’s  lemma  tells  us  that  e®  is  a  d-edge,  t.e.,  the  lines  containing  e^_x  and 
e^+1  intersect  on  P’s  side  of  the  line  containing  ek.  This  implies  (ot+2  —  a*_i)  mod  2n  >  ir, 
and  this  inequality  can  hold  for  at  most  three  values  of  k. 

To  identify  those  intervals  that  might  contain  the  contact  point  of  the  non-flush  edge  of 

Q,  we  need  only  compute  ak  for  each  interval  lk  and  then  (ak+2  —  ofc_i)  mod  2 ir  for  each  k, 
which  may  be  done  in  0(d)  time.  ■ 

We  can  check  the  possibility  of  the  nonflush  edge  lying  in  the  interval  Ik  as  follows.  Let 
ef  be  any  edge  of  P  in  Ik-\-  Let  ef  be  any  edge  of  P  in  J*+J.  Define  QtiJ  to  be  the  set  of 
all  circumscribing  d- gons  Q  satisfying  the  following  constraints: 

1.  Q's  (k  —  l)-st  contact  point  is  the  edge  ef , 

2.  Q' s  fc-th  contact  point  is  an  edge  or  vertex  in  Ik, 

3.  Q's  fc-th  edge  is  a  d-edge, 

4.  Q's  (k  +  l)-st  contact  point  is  the  edge  ef ,  and 

5.  for  1  <  l  <  k  —  1  and  k  +  1  <  /  <  d,  Q's  /-th  contact  point  is  an  edge  in  /(. 

Also  define 

a  =  Ua.,- 

i 

Lemma  3.7  For  any  ef  €  Ik-i,  we  can  find  a  circumscribing  d-gon  Q,  €  Q,  of  minimal 
area  in  0(n)  time. 


Proof  We  begin  with  a  few  more  definitions  (borrowed  from  [ACY85]).  An  h-sided  (i,  j)- 
chain  is  a  polygonal  chain  C  =  {ef . . . . ,  ef}  such  that 

1.  ef  is  flush  with  ef, 

2.  ef  is  flush  with  ef ,  and 

3.  for  1  <  /  <  h ,  ef  has  a  contact  point  in  7(*+/)modd- 

C  is  flush  if  all  its  edges  are  flush.  The  extra  area  of  C  is  the  area  of  the  bounded  (but 
perhaps  disconnected)  region  between  C  and  P. 

For  each  edge  ef  £  h+i,  let  Cfj  be  an  optimal  flush  (d—  l)-sided  (j,  i)-chain  ( i.e its  extra 
area  is  minimal),  and  let  C-j  be  an  optimal  3-sided  (*,  j)-chain.  Combining  Ctj  and  C-  j  gives 
am  optimal  circumscribing  d- gon  Qij  from  Qjj.  Moreover,  given  Qitj  for  each  ef  £  h+ 1,  we 
can  pick  the  best  of  these  d-gons  in  0(n)  time  to  obtain  an  optimal  circumscribing  d-gon 
from  Qi.  Thus,  if  we  can  find  an  optimal  flush  (d  —  l)-sided  (,;,?)-chain  and  an  optimal 
3-sided  chain  for  each  ef  £  Ik+X  in  O(n)  time,  we  will  have  established  the  lemma. 

To  find  an  optimal  flush  (d  —  l)-sided  (j,  i)  chain  for  each  ef  £  J*+i,  we  first  define  the 
(d  —  2)-dimensional  path-decomposable  array  A  =  . ,d_2},  where 

=  *^»1.«2  "h  "b  *  ‘  '  w*d-k-\.id-k  "b  W‘d-k+l,‘d-k  W'd-k+ 1  .<<-*+2  "b  ’  '  '  *b  ^id_2,f 

(We  defined  the  arrays  W  =  and  W  =  earlier  in  this  section,  in  the  context 

of  optimal  flush  circumscribing  d-gons.)  Now  suppose  7/  =  \jhji+\]  for  1  <  /  <  d.  Let 
m  =  ji+i  -  ji  +  1,  and  let  A'  be  the  nx  x  n2  x  -  •  •  x  nd-2  subarray  of  A  containing  entries 

a„,i2 . id_3  that  satisfy  j{i+k)m°dd  <  U  <  J(i+fc+i)modd  for  1  <  /  <  d  -  2.  The  optimal  flush 

(d  —  l)-sided  (j,  i)  chain  corresponds  to  the  maximum  entry  in  the  j-th  plane  of  A!  (the 
plane  containing  those  entries  of  A!  whose  first  coordinate  is  j).  Since  JZf",2  n,  =  O(n), 
Theorem  2.10  tells  us  that  we  can  compute  these  plane  maxima  in  O(n)  time. 

To  find  an  optimal  3-sided  (i,  j)-chain  for  each  ef  €  Ik+i,  we  first  recall  that  Lemma  3.4 
tells  us  we  need  only  consider  3-sided  chains  whose  middle  edge  is  a  balanced  d-edge.  We 
will  call  this  middle  edge  a  closing  edge  for  ef  and  ef . 

Now  suppose  ef  is  a  balanced  closing  d-edge  for  ef  €  Ik-i  and  ef  €  W  e?’s  first 
endpoint  must  lie  on  the  line  Li  containing  ef ,  and  its  second  endpoint  must  lie  on  the 
line  Lj  containing  ef.  Moreover,  since  ef  is  balanced,  its  second  endpoint  must  also  lie 
on  the  chain  C  defined  as  follows.  For  each  ef  €  /*,  let  efj.i  be  the  segment  on  the  line 
containing  ef  whose  first  and  second  endpoints  are  twice  as  far  from  X,-  as  are  ef ’s  first  and 
second  endpoints,  respectively.  C  consists  of  these  segments,  plus  segments  eft  parallel  to  i, 
connecting  eft_l  and  ef1+i  for  1  <  /  <  \Ik\.  This  is  suggested  in  Figure  3.4.  We  will  denote 
the  endpoints  of  ef  by  vf  and  vf+1 . 

Using  what  we  know  about  the  angles  of  rays  containing  ef  €  C  and  ef  £  Ik+X  form 
with  respect  r'cw,  we  can  prove  the  following  two  propositions  about  C  and  the  lines  Lj 
containing  edges  in  w 

Proposition  3.8  For  all  ef  £  Ik+t,  Lj  intersects  C  at  most  once. 

Proposition  3.9  For  all  ef  and  ef  in  Ik+i,  j'  >  j,  Lj  intersects  C  before  Lj>  does,  i.e.,  if 
Lj  intersects  ef ,  then  Lj<  intersects  ef, ,  l'  >  l. 
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Figure  3.4:  Every  balanced  closing  d-edge  for  a  particular  ef  G  h- 1  and  any  ef  G  h+\  must 
have  one  endpoint  on  the  chain  C. 

The  first  proposition  tells  us  that  there  is  at  most  one  balanced  closing  d-edge  corre¬ 
sponding  to  each  e f  G  h+i-  The  second  allows  us  to  find  all  these  closing  edges  in  0(n) 
time.  We  begin  by  computing  C,  which  can  be  done  in  0(n )  time.  We  then  check  whether 
the  line  L  containing  the  first  edge  in  /k+i  intersects  ef .  If  L  does  not  intersect  this  segment, 
we  check  whether  it  intersects  ef .  If  again  there  is  no  intersection,  we  move  on  to  ef .  We 
continue  in  this  manner  until  we  find  the  segment  ef  that  L  intersects.  This  gives  us  the 
closing  edge  for  the  first  edge  of  Jfc+1  Next,  we  check  whether  the  line  V  containing  the 
second  edge  of  Ik+\  intersects  ef.  If  there  is  no  intersection,  we  move  on  to  segment  ef+v 
We  continue  in  this  manner  until  all  closing  edges  are  found.  It  is  easy  to  see  that  only 
0(n )  time  is  required.  Also,  Proposition  3.9  guarantees  that  this  approach  will  find  all  of 
the  desired  intersections,  which  in  turn  gives  us  the  optimal  3-sided  (t,,;)-chains.  ■ 

To  complete  our  algorithm,  we  require  a  third  corollary  to  Theorem  3.2. 

Corollary  3.10  For  any  ef  and  ef  in  Ik- 1,  every  optimal  Q,  €  Q,  interleaves  every  optimal 

Qi'  €  Qi>. 

Proof  This  corollary  follows  immediately  from  Theorem  3.2,  since  Q,  and  are  closed 
under  edge  exchange.  ■ 

This  corollary,  together  with  Lemma  3.7,  allows  us  to  find  an  optimal  circumscribing 
d-gons  for  each  Qi  such  that  ef  €  h- 1  in  O(nlgn)  total  time.  We  use  the  natural  divide- 
and-conquer  approach  of  Theorem  2.6.  Let  /;  =  \ji,ji+\]  for  1  <  /  <  d,  and  let  = 
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jt+ 1  —  ji  -f  1.  We  first  find  an  optimal  circumscribing  d-gon  Qi  €  Qi  for  i  =  +  fn*_i/2"|. 

By  Corollary  3.10,  the  contact  point  of  Qi  in  interval  //,  1  <  /  <  d,  splits  that  interval  into 
two  intervals  7/  and  7",  such  that  the  contact  points  of  any  optimal  circumscribing  d- gon  for 
£>,•»,  jk- i  <  <  i,  must  lie  in  the  intervals  7j, . . . ,  7j,  and  the  contact  points  of  any  optimal 

circumscribing  d- gon  for  £,■»,  i  <  i'  <  jk ,  must  lie  in  the  intervals  7", . . . ,  I'j.  If  we  recursively 
solve  the  two  subproblems  associated  with  the  intervals  7{ , . . . ,  I'd  and  the  intervals  I", ... ,  I'j, 
we  obtain  a  recurrence  with  solution  0(n  lgn^)  =  O(nlgn)  for  the  time  required  to  find 
an  optimal  circumscribing  d-gon  for  each  Q,  such  that  ef  €  Jk-i- 

By  choosing  the  best  of  the  n^-i  circumscribing  d-gons  obtained  in  this  manner  (which 
can  be  done  in  0(n)  time),  we  obtain  a  minimum  area  circumscribing  d-gon;  thus, 

Theorem  3.11  Given  a  convex  n-gon,  a  minimum  area  circumscribing  d-gon  can  be  com¬ 
puted  in  0(dn  +  n  lg  n)  time. 

3.2  Perimeter  Minimization 

Just  as  the  techniques  given  by  Boyce  et  al.  [BDDGS5]  for  finding  a  maximum  perimeter 
inscribed  d-gon  can  be  used  to  find  a  minimum  area  flush  circumscribing  d-gon  in  0(dn  + 
nlgn)  time,  they  can  also  be  used  to  find  a  minimum  perimeter  flush  circumscribing  d-gon 
in  the  same  amount  of  time.  Unfortunately,  we  lack  a  lemma  analogous  to  Lemma  3.4  for  the 
case  of  perimeter  minimization,  and  thus  it  is  unclear  how  to  use  this  minimum  perimeter 
flush  circumscribing  d-gon  to  obtain  a  minimum  perimeter  arbitrary  circumscribing  d-gon.  In 
fact,  the  only  previous  algorithm  for  the  minimum  perimeter  circumscribing  d-gon  problem, 
due  to  DePano  [DeP87],  is  restricted  to  the  case  of  d  =  3.  Specifically,  DePano  obtained  an 
<9(n3)  time  algorithm  for  the  minimum  perimeter  circumscribing  triangle  problem.  In  doing 
so,  he  proved  the  following  two  lemmas,  which  we  will  use  in  obtaining  an  O(nlgrc)  time 
algorithm  for  the  problem. 

Lemma  3.12  ([DeP87])  For  every  convex  n-gon  P,  there  exists  a  minimum  perimeter 
circumscribing  triangle  Q  with  at  least  one  edge  flush  with  P. 

Lemma  3.13  ([DeP87])  For  any  triple  ( i,j ,  k),  there  is  a  unique  minimum  perimeter  cir¬ 
cumscribing  triangle  Tij,k  whose  first  edge  ef  is  flush  with  edge  ef  of  P,  whose  second  edge 
ef  contains  vertex  vf  of  P,  and  whose  third  edge  ef  contains  vertex  vf  of  P.  Moreover,  this 
triangle  satisfies  the  following  condition:  there  exists  a  point  p  on  ef  and  P  and  a  point  q 
on  ef  and  P  such  that  the  distance  between  p  and  the  endpoint  shared  by  ef  and  ef  equals 
the  distance  between  q  and  the  endpoint  shared  by  ef  and  ef  (see  Figure  S.5).  Note  that  p 
must  be  vf  if  ef  is  not  flush  with  P,  and  q  must  be  vf  if  ef  is  not  flush  with  P. 

Lemma  3.12  tells  us  that  we  need  only  consider  circumscribing  triangles  from  the  set 

T  =  {7U*  |  1  <  *  <  n,  1  <  j  <  k  <  n}, 

as  this  set  must  contain  a  minimum  perimeter  circumscribing  triangle.  For  any  t,  j ,  and  k, 
Lemma  3.13  allows  us  to  compute  Tjj,*  (and  its  perimeter)  in  constant  time  (see  [DeP87]  for 
details). 
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Figure  3.5:  If  Aabc  is  the  minimum  perimeter  circumscribing  triangle  containing  edge  ef  and 
vertices  vf*  and  vk  of  P,  then  there  exist  points  p  and  q  on  both  P  and  Aabc  such  that  the  distance 
li  between  a  and  p  equals  the  distance  I2  between  b  and  q. 

DePano’s  0(n3)  time  algorithm  for  the  minimum  circumscribing  triangle  problem  is  noth¬ 
ing  more  than  a  brute  force  search  through  T  —  he  just  computes  the  perimeter  of  every 
triangle  in  T  and  then  outputs  one  with  minimal  perimeter.  We  obtain  a  substantially  better 
algorithm  by  identifying  a  multidimensional  Monge  array  associated  with  T  and  applying 
our  array-searching  techniques. 

Specifically,  we  define  the  following  n  x  (2n  — 2)  x  (2n  — 1)  array  A  =  {al]k}.  For  1  <  1  <  n 
and  i  <  j  <  k  <  i  +  n,  we  let  aljk  be  the  perimeter  of  Ttj^  (which  we  denote  per(7’i>Jit)), 
provided  this  triangle  exists  (the  indices  j  and  k  are  modulo  n).  If  this  triangle  does  not 
exist,  we  let  ai}k  be  00.  For  all  other  values  of  t,  j,  and  k ,  we  also  define  a, }k  to  be  00. 

Now,  every  triangle  in  T  is  represented  in  A,  and  the  only  entries  in  A  that  are  less  than 
00  correspond  to  triangles.  Thus,  by  Lemma  3.12,  the  minimum  entry  in  A  corresponds  to 
the  minimum  perimeter  circumscribing  triangle.  Moreover,  since  Lemma  3.13  implies  each 
entry  of  A  can  be  computed  in  constant  time,  the  following  lemma  gives  us  an  0{n  Ig  n)  time 
algorithm  for  computing  a  minimum  perimeter  circumscribing  triangle. 

Lemma  3.14  A  satisfies  the  Monge  condition. 

Proof  Consider  any  two-dimensional  plane  of  A  corresponding  to  a  fixed  value  i  of  yl’s 
first  coordinate.  Suppose  this  plane  does  not  satisfy  the  Monge  condition,  t.e.,  there  exist 
j  <  j'  and  k  <  k'  such  that 

>  aijk'  +  O’ij'k-  (3-1) 

We  cannot  have  j'  >  k,  since  this  would  mean  a,y*  =  00,  which  contradicts  (3.1).  Similarly, 
Tijk>  must  exist,  which  implies  T,y*,  T,yk>,  and  Tiyk  also  exist.  Now  suppose  we  exchange 
edges  between  Tijk>  and  TtJ<k.  Specifically,  we  let  Vi*  be  the  triangle  formed  by  the  first  and 
second  edges  of  T,jk.  and  the  third  edge  of  T,yk  (extended  or  shortened,  as  necessary,  to  form 
a  triangle),  and  let  T{yk,  be  the  triangle  formed  by  the  first  and  second  edges  of  X,y*  and 
the  third  edge  of  This  is  suggested  in  Figure  3.6.  Now  it  is  easy  to  verify  that 

per(Tijfc)  +  per(7’yfc.)  <  per(7,iifc')  +  P^Xyt). 


22 


Figure  3.6:  Suppose  we  exchange  edges  between  the  triangles  and  Tjy*  (indicated  by 

the  dotted  lines)  and  obtain  the  triangles  and  T'j,k,  (indicated  by  the  shaded  lines).  Then 
per(T/jfc)  -f  per(T/yfc,)  <  per(T,J*.-)  -f  per(T,y«t),  since  the  length  of  the  line  segment  a  plus  the 
length  of  the  line  segment  6  is  at  least  the  length  of  the  line  segment  b  plus  the  length  of  the  line 
segment  c. 

Furthermore,  permit)  <  per(7^jfc)  and  per(7i;<*<)  <  per(T'yk> ),  since  T,jk  and  T.yk1  have 
minimal  perimeter  among  circumscribing  triangles  that  share  their  contact  points.  Thus, 

per(r.jfc)  +  per(T,y*')  <  per(7V,*<)  +  per(r,7-*), 
which  contradicts  (3.1). 

In  an  similar  fashion,  we  can  show  that  every  two-dimensional  plane  of  A  corresponding 
to  a  fixed  value  of  the  second  coordinate  or  a  fixed  value  of  the  third  coordinate  also  satisfies 
the  Monge  condition.  ■ 

It  is  unclear  whether  A  is  Monge-composite,  since  fixing  any  two  indices  of  A  does  not 
fix  the  both  the  edges  of  the  circumscribing  triangle  corresponding  to  these  indices.  Never¬ 
theless,  Lemma  3.14  allows  us  to  apply  the  array-searching  algorithm  given  in  Theorem  2.6 
and  obtain: 

Theorem  3.15  Given  a  convex  n-gon  P,  a  minimum  perimeter  triangle  Q  circumscribing 
P  can  be  computed  n  O(nlgn)  time. 

It  remains  open  whether  a  minimum  perimeter  circumscribing  triangle  can  be  found  in 
o(n  lgn)  time  and  whether  there  exist  efficient  algorithms  for  computing  minimum  perimeter 
circumscribing  d- gons. 


4  Dynamic  Programming  Using  Monge  Conditions 

4.1  Frances  Yao’s  Paper  Revisited 

In  [YaoSO],  Yao  considered  the  following  dynamic  programming  problem.  Let  w(i,j)  be  a 
weight  function  defined  for  1  <  i  <  j  <  n  and  satisfying  the  following  two  constraints:  for 
1  <  i  <  i'  <  j  <  j'  <  n, 


w{i,j)  +  u>(i\  j')  <  w(i,j')  +  (4.1) 

and 

w(i',j)  <  w(i,j').  (4.2) 

(Yao  called  the  first  constraint  the  quadrangle  inequality  and  referred  to  functions  satisfying 
the  second  constraint  as  monotone  on  the  lattice  of  intervals,  ordered  by  inclusion.)  We  want 
to  compute  c(t,j)  for  1  <  i  <  j  <  n,  where  c(i,  i)  =  0  for  1  <  i  <  n  and 

c[i,j)  =  tv(ij)  +  min  {c{i.k  -  1)  +  c(k,j)} 

t<K; 

for  1  <  i  <  j  <  n.  In  [YaoSO],  Yao  showed  that  this  problem  can  be  solved  in  0(n2)  time.  She 
used  this  result  to  provide  a  simple  alternate  solution  to  Knuth’s  problem  regarding  optimal 
binary  trees  [I\nu73]:  she  also  extended  this  approach  to  the  computation  of  optimal  t-ary 
trees. 

In  this  subsection,  we  reformulate  Yao’s  dynamic  programming  problem  in  terms  of 
multidimensional  Monge  arrays  and  then  use  a  result  due  to  Wilber  [WilSS],  to  obtain  an 
alternate  0{n2)  time  solution  for  the  problem.  We  begin  with  a  lemma  proved  by  Yao  in 
obtaining  her  0{n2)  time  bound. 

Lemma  4.1  ([Yao80])  The  cost  function  c(i,j)  satisfies  the  quadrangle  inequality,  i.e., 
c(i,j)  +  c(i',j')  <  c(ij')  +  c(i',j)  for  1  <  i  <  i’  <  j  <  j '  <  n. 


Yao’s  quadrangle  inequality  is  precisely  the  Monge  condition,  except  that  the  functions 
w(i,j)  and  c(i,j)  do  not  correspond  to  complete  arrays.  However,  if  we  let  W  — 
denote  the  n  x  n  array  where 


f  u>(i,  j)  if  i  <  j , 

(  oo  otherwise, 


and  we  let  C  =  {c,j}  denote  the  n  x  n  array  where 


f  c(t,j)  if  i  <  j, 

(  oo  otherwise, 


then  both  W  and  C  satisfy  the  Monge  condition.  (Note  that  it  is  important  that  we  define 
Wij  and  Cij  to  be  oo  when  i  >  j  —  if  we  instead  define  Wij  and  c,  j  to  be  — oo  when 
i  >  j,  as  we  might  want  to  do  if  we  were  maximizing  instead  of  minimizing,  then  W  and 
C  would  not  satisfy  the  Monge  condition  or  the  inverse  Monge  condition.)  Yao’s  dynamic 
programming  problem  then  boils  down  to  computing  the  entries  of  C .  Unfortunately,  the 
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results  of  [AKM*S7]  seems  inapplicable  at  this  point,  at  least  in  a  straightforward  manner, 
as  we  are  neither  interested  in  the  row  minima  of  C  nor  are  the  entries  of  C  readily  available. 

Now  consider  the  three-dimensional  array  D  =  where  we  define  d,tk,3  =  + 

Cfc  j  +  Wjj.  Since  we  now  have  c(jJ  =  minj</i<ri  we  have  transformed  the  problem  of 
computing  the  entries  of  C  to  the  problem  of  computing  the  tube  minima  of  D.  Moreover,  D 
is  clearly  a  path-decomposable  Monge-composite  array  (which  means  it  satisfies  the  Monge 
condition),  but  again  the  array’s  entries  are  not  directly  available.  At  this  point,  however, 
we  can  apply  the  results  of  Wilber  [WilSS]. 

Wilber  considered  the  concave  least-weight  subsequence  problem:  given  a  weight  function 
v(p,  q)  defined  for  0  <  p  <  q  <  m  and  satisfying  what  Yao  calls  the  quadrangle  inequality, 
compute  f(q)  for  0  <  q  <  m,  where  /( 0)  =  0  and  f(q)  =  mino<p<9{/(p)  -f  v(p,  q)}  for 
l  <  q  <  m.  We  can  convert  this  problem  to  an  array- searching  problem  by  letting  V  =  {uPi,} 
denote  the  (m  +  1)  x  (m  +  1)  array  where 


■w  -  { 


q)  if  P  <  9< 
otherwise, 


and  letting  E  =  {ePi,}  denote  the  (m  +  1)  x  (m-t-1)  array  where  epq  =  f(p)+vPil).  The  concave 
least-weight  subsequence  problem  is  now  the  column  minima  problem  for  E.  Moreover,  since 
v(p,q)  satisfies  the  quadrangle  inequality,  V  satisfies  the  Monge  condition,  which  implies  E 
satisfies  the  Monge  condition.  In  [WilSS],  Wilber  showed  that  this  problem  can  be  solved  in 
0(n)  time,  even  though  computing  an  arbitrary  entry  of  E  is  not  a  constant  time  operation. 

Returning  to  Yao’s  dynamic  programming  problem,  consider  the  two-dimensional  plane 
of  D  corresponding  to  a  particular  value  5  of  the  first  coordinate.  The  tube  minima  of  D 
lying  in  this  plane  are  just  the  column  minima  of  this  plane.  Moreover,  the  column  minima 
problem  for  the  plane  is  an  instance  of  the  concave  least-weight  subsequence  problem.  In 
particular,  m  =  n  —  i,  vP),  =  Cp+i+1,,+1  +  u>il9+i,  and  f(q)  =  Cj,g+ Note  that  the  array  V 
satisfies  the  Monge  condition  because  the  array  C  satisfies  the  Monge  condition.  Provided 
ctj  is  known  for  i  >  i  (so  that  vp<q  can  always  be  computed  in  constant  time),  we  can  solve 
this  problem  in  O(m)  time  using  the  algorithm  given  by  Wilber  in  [Wil8fh  Thus,  to  compute 
the  tube  minima  of  Z),  we  need  only  apply  Wilber’s  algorithm  n  time:,  hi.  >  the  plane 
corresponding  t  o?  =  n,  then  to  the  plane  corresponding  to  i  =  n  —  1,  and  so  on  down  to  the 
plane  corresponding  to  i  —  1.  This  gives  the  following  theorem. 

Theorem  4.2  If  the  weight  function  w{i,j)  satisfies  (4.1)  and  (4-2),  then  Yao’s  dynamic 
programming  problem  can  be  solved  in  0(n2)  time. 

Now  suppose  the  weight  function  w(i,j)  satisfies 

w(ij)  +  w(i',j')  >  w(i,j')  +  w(i\j)  (4.3) 


w(i',j)  >  iv(ij')  (4.4) 

for  1  <  i  <  i'  <  j  <  f  <  n.  By  adopting  the  procedure  given  above,  we  can  obtain  the 
following  theorem. 


Theorem  4.3  If  the  weight  function  w(i,j)  satisfies  (4-3)  and  (4-4)-  then  Yao's  dynamic 
programming  problem  can  be  solved  in  0(n2a(n))  time,  where  a( n)  denotes  the  inverse  Ack- 
ermann  function. 

Proof  This  proof  is  identical  to  that  given  above  except  that  now  the  arrays  IV,  C,  and 
D  defined  above  are  no  longer  Monge  (or  even  totally  monotone).  The  column  minima 
problem  for  the  plane  of  D  corresponding  to  a  particular  value  i  of  the  first  coordinate 
is  now  equivalent  to  the  convex  least-weight  subsequence  problem  considered  in  [EGG88]: 
given  a  weight  function  v(p,q )  defined  for  0  <  p  <  q  <  m  such  that  v{p,  q)  4-  v{p\  <f)  > 
v(p,q’)  +  v(p',q)  for  0  <  p  <  p'  <  q  <  q'  <  m,  compute  f(q)  for  0  <  q  <  m,  where  /(0)  =  0 
and  f  (q)  =  min0<p<g{/(p)  +  u(p,  q)}  for  1  <  q  <  m.  Klawe  and  Kleitman  [KK88]  have  shown 
that  this  problem  can  be  solved  in  0(na(n))  time.  Thus,  by  applying  Klawe  and  Kleitman ’s 
algorithm  n  times,  we  can  obtain  the  entries  of  C  in  0{n2a(n))  time.  ■ 

It  remains  open  whether  the  time  complexity  given  in  Theorem  4.3  can  be  improved  from 
0(n2a(n))  to  0(n2). 

4.2  Waterman’s  Problem 

The  primary  structure  of  a  single-stranded  RNA  molecule  is  the  sequence  of  ribonucleotides 
or  bases  making  up  the  molecule.  When  the  primary  structure  of  an  RNA  molecule  is 
known,  the  question  of  which  bases  form  pairs  (thus  pinching  off  loops  in  the  chain  of  bases) 
becomes  important  —  this  structure  is  referred  to  as  the  secondary  structure  of  the  RNA. 
In  [Wat7S],  Waterman  argued  that  predicting  the  secondary  structure  of  an  RNA  molecule 
from  its  primary  structure  is  closely  related  to  solving  some  dynamic  programming  problems. 
Among  the  dynamic  programming  problems  he  described  is  the  following:  given  a  weight 
function  w(i,j)  defined  for  1  <  i  <  j  <  n  and  satisfying  the  quadrangle  inequality  (re., 

w(k,  l)  +  w{p,  q)  <  w(k ,  q)  +  w(p ,  /)  (4.5) 

for  1  <  k  <  p  <  /  <  q  <  n),  compute 

e(*,j)  =  nun.  -f  w(i'  -f  j',i  +  j) 

i<r<t 

for  all  i  and  j  between  1  and  n,  where  e(t,  1)  and  e(l,j)  are  given  for  all  i  and  j  and  c(i,j) 
can  be  computed  from  e(t ,j)  in  constant  time. 

The  naive  approach  gives  an  0(n4)  time  algorithm  for  this  dynamic  programming  prob¬ 
lem.  Waterman  and  Smith  improved  this  time  bound  to  0(n3)  in  [WS86],  and  then  recently 
Eppstein,  Galil,  and  Giancarlo  [EGG88]  obtained  an  0(n2  lg2  n)  time  algorithm  for  the 
problem.  (They  also  obtained  an  0(n2  lg  nig  lg  n)  time  algorithm  for  the  special  case  of  w  a 
logarithm  or  other  simple  function.)  Using  array-searching,  we  show: 

Theorem  4.4  If  the  weight  function  w(i,j)  satisfies  the  quadrangle  inequality  (4-5),  then 
Waterman’s  dynamic  programming  problem  can  be  solved  in  0(n2lgn)  time. 
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Proof  We  only  prove  the  theorem  for  the  case  c(i,j)  =  e(i,j ),  but  the  proof  is  easily 
generalized  to  any  c(i,j)  that  can  be  computed  from  e(i,j)  in  constant  time. 

We  begin  by  extending  the  function  w(i,j)  to  form  an  n  x  n  array,  just  as  we  did  in  the 
context  of  Yao’s  dynamic  programming  problem.  In  particular,  we  let  W  =  {u?,-^}  denote 
the  array  where 

w{i,j)  if  i<j, 
oc  otherwise. 


U>; 


-{ 


As  before,  this  array  satisfies  the  Monge  condition.  We  will  also  view  the  function  e(i,  j)  as 
an  n  x  n  array  E  =  {e,j}  where  a,j  =  e(i,  j),  so  that  the  problem  is  now  that  of  computing 
the  entries  of  E. 

We  now  require  a  few  definitions  and  observations.  We  will  say  that  an  entry  e,j  from  E 
strictly  dominates  an  entry  e{>^  if  i  >  i'  and  j  >  j'.  Note  that  an  entry  of  E  depends  only 
on  those  entries  it  strictly  dominates.  We  will  also  define  the  k-th  diagonal  of  E  to  consist 
of  those  entries  e,j  such  that  ?  +  j  =  k.  Note  that  the  only  entry  on  the  k-th  diagonal  of  E 
that  we  need  to  consider  in  computing  e,0  is  the  minimum  entry  e, on  the  diagonal  that 
is  strictly  dominated  by  et>J.  We  will  refer  to  this  entry  as  a  diagonal  minimum. 

With  this  is  mind,  we  define 


dk(*,j)  =  min  e,v. 

i  <i 

j'<i 

i'+j'-k 

(If  k  =  1  or  k  >  i  +  j  —  2,  we  let  d*(i,  j)  =  oo.)  We  then  have 

e.,j  =  mmjdk(tj)  +  wk^}. 

We  can  now  describe  our  algorithm  for  computing  the  entries  of  E.  We  use  a  divide- 
and-conquer  approach,  reminiscent  of  an  algorithm  given  by  Aggarwal  and  Suri  [ASST]  for 
finding  the  largest  empty  corner  rectangle.  Without  loss  of  generality,  we  assume  n  is  a 
power  of  2.  We  begin  by  partitioning  E  into  four  n/2  x  n/2  subarrays  A,  B,  C ,  and  D. 
A  contains  entries  from  the  first  n/2  rows  and  first  n/2  columns  of  E ,  B  contains  entries 
from  the  first  n/2  rows  and  last  n/2  columns,  C  contains  entries  form  the  last  n/2  rows 
and  first  n/2  columns,  and  D  contains  the  remaining  entries.  We  then  recursively  compute 
the  entries  of  A;  this  requires  T(n/ 2)  time,  where  T(n)  is  the  time  required  to  compute  the 
entries  of  an  n  x  n  array.  Next,  we  compute  the  effect  of  A  on  B,  i.e.,  for  each  entry  of  J9, 
we  obtain  an  upper  bound  on  its  value  based  on  the  entries  of  A.  (Each  entry  in  B  is  the 
minimum  of  a  number  of  terms;  some  of  these  terms  depend  only  on  entries  from  B ;  the 
rest  depend  only  on  entries  from  A;  we  compute  the  minimum  of  those  terms  dependent  on 
entries  from  A.)  We  will  explain  later  how  this  may  be  accomplished  in  0(n2)  time.  We 
then  recursively  compute  the  entries  of  B  in  T(n/2 )  time,  taking  into  account  the  effect 
of  A  on  B,  i.e .,  in  assigning  an  entry  of  B  its  value,  we  always  take  the  minimum  of  the 
value  we  have  recursively  computed  for  this  entry  based  on  the  entries  of  B,  and  the  value 
we  have  computed  for  this  entry  based  on  the  entries  of  A.  We  repeat  this  process  for  C, 
computing  the  effect  of  A  on  C  in  0(n2)  time  and  then  recursively  computing  C  in  T(n/ 2) 
time.  Finally,  we  compute  the  effect  of  A  on  D ,  the  effect  of  B  on  D ,  and  the  effect  of  C  on 
Z?,  all  in  0(n2)  time,  and  then  recursively  compute  D  in  T(nj 2)  time. 
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This  yields  the  recurrence  T(n)  <  4!T(n/2)  -(-  0(n2)  for  the  time  to  compute  all  entries 
of  an  n  x  n  array.  Since  T(l)  =  0(1),  T(n)  =  0(n2lgn). 

We  can  compute  the  effect  of  A  on  B  as  follows  (computing  the  effect  of  A  on  C ,  A  on 
D ,  B  on  D ,  and  C  on  D  may  be  done  in  an  analogous  manner).  For  1  <  i  <  n/2,  each  entry 
ei,j+n/ 2  in  row  i  of  B  strictly  dominates  the  same  set  of  entries  in  A  and  thus  depends  on 
the  same  diagonal  minima.  This  motivates  defining  the  (i  +  n/2)  x  n/2  array  X'  = 
where  =  d*(i,  1  +  n/2)  +  Wk,i+j+n/ 2-  e«,j+n/2  can  then  be  expressed  as  the  minimum  of 
the  following  two  minima: 


mm 

l<A<t+n/2 


XkJ 


and 

min  {e;<j»+n/2  +  u),<+j-+„/2,<+j+n/2}- 

i<«  <« 

As  the  first  of  these  quantities  depends  only  on  entries  in  A  and  the  second  depends  only  on 
entries  in  B ,  computing  the  effect  of  A  on  B  reduces  to  computing  the  column  minima  in 
all  n  of  the  X'.  Furthermore,  observe  that  for  all  i  between  1  and  n/ 2,  the  array  Xx  satisfies 
the  Monge  condition,  since  W  satisfies  the  Monge  condition.  Thus,  we  have  the  following 
lemma. 

Lemma  4.5  The  column  minima  of  X'  for  all  i  between  1  and  n/2  can  be  computed  in 
0(n2)  time. 


Proof  of  Lemma  4.5  The  <4(i,n/2)  for  1  <  i  <  n/2  and  1  <  k  <  n  can  be  computed  in 
0(n2)  time,  since  <4(1,  n/2)  =  oo  for  1  <  k  <  n, 


f  min{dfc(i  -  l,n/2),  if  i  <  k  <  n/2  +  i  -  2, 

|  4(i-l,n/2)  otherwise, 


for  2  <  i  <  n/2  and  1  <  k  <  n,  and  eitj  is  known  for  i  and  j  between  1  and  n/2.  Once  this 
has  been  done,  any  x'k<:  can  be  computed  in  constant  time,  which  means  we  can  apply  the 
algorithm  of  [AKM*87]  n  times  and  obtain  all  the  desired  column  minima  in  an  additional 
0(n2)  time.  ■ 


This  completes  our  proof  of  Theorem  4.4.  ■ 


In  [EGG88],  Eppstein,  Galil,  and  Giancarlo  also  considered  the  variant  of  Waterman’s 
problem  in  which  the  weight  function  w(i,j)  satisfies  the  inverse  quadrangle  inequality,  t.e., 

w(k,  l)  +  w(p,  q)  >  w(k ,  g)  +  w(p,  l )  (4.6) 

for  Jfc  <  p  <  /  <  q).  They  obtained  an  0(n2lg2n)  time  algorithm  for  this  problem  using 
the  same  techniques  they  employed  for  weight  functions  satisfying  the  quadrangle  inequality. 
We  can  apply  our  techniques  to  this  variant  of  Waterman’s  problem,  too,  and  obtain  the 
following  theorem. 


Theorem  4.6  If  the  weight  function  w(i,j)  satisfies  the  inverse  quadrangle  inequality  (4-6), 
then  Waterman’s  dynamic  programming  problem  can  be  solved  in  0(n2lgn)  time. 
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Proof  This  proof  is  identical  to  that  given  above  for  weight  functions  satisfying  the  quad¬ 
rangle  inequality,  except  that  the  array  W  defined  above  is  no  longer  Monge  (or  even  totally 
monotone).  However,  the  arrays  X'  are  Monge  (they  satisfy  the  inverse  Monge  condition), 
since  they  are  formed  from  subarrays  of  Hr  that  are  Monge.  Thus,  we  can  apply  the  algorithm 
of  [Al\M*87j  as  before.  ■ 

Note  that  the  divide-and-conquer  approach  our  algorithms  employ  can  be  used  to  obtain 
an  O(nlgn)  time  solution  to  the  one-dimensional  version  of  Waterman’s  problem,  namely, 
the  concave  or  convex  least-weight  subsequence  problem  considered  by  [HL87,  KK88,  Wil88]. 
However,  Wilber  [Wil88]  solves  the  concave  least- weight  subsequence  problem  in  0(n)  time 
by  an  elegant  modification  of  [AKM*87],  and  Klawe  and  Kleitman  [KK88]  solve  the  convex 
least-weight  subsequence  problem  in  0(na(n))  time.  This  suggests  that  perhaps  an  o(n2  lg  n) 
solution  to  Waterman's  problem  might  be  obtainable. 


5  VLSI  River  Routing 

In  the  section,  we  apply  array-searching  to  a  number  of  problems  that  arise  in  VLSI  river 
routing.  These  problems  involve  two  pieces  (or  modules)  P  and  Q  of  a  VLSI  circuit  that 
need  to  be  wired  together.  We  model  P  as  a  sequence  of  n  points  (or  terminals)  lying  on  a 
horizontal  line  in  the  plane.  We  let  p i  denote  both  the  ?‘-th  terminal  of  P  and  the  terminal’s 
distance  from  px.  We  also  assume  pi  <  P2  <  •  ■  ■  <  Pn-  Similarly,  we  model  Q  as  another 
sequence  of  n  terminals  lying  on  a  second  horizontal  line,  where  <7,  denotes  both  the  i-th 
terminal  of  Q  and  the  terminal’s  distance  from  ft,  and  ft  <  g2  <•••<?„•  A  placement 
for  Q  (relative  to  P)  consists  of  an  offset  and  a  separation ,  which  give  ft’s  horizontal  and 
vertical  position,  respectively,  relative  to  pt.  Given  a  placement  for  Q,  a  routing  is  a  set  of 
n  wires  (represented  by  continuous  curves),  where  the  i-th  wire  connects  p,  to  qt.  In  river 
routing  (as  opposed  to  more  general  wiring  models),  these  wires  must  be  nonintersecting. 
Furthermore,  the  wires  are  constrained  to  satisfy  a  set  of  technology-dependent  design  rules. 
At  a  minimum,  the  distance  between  any  two  wires  is  always  at  least  some  fixed  constant, 
which  we  take  to  be  1.  Examples  of  other  constraints  include  requiring  that  the  wires  are 
rectilinear,  i.e.,  they  lie  on  a  rectilinear  integer  grid  (as  they  do  in  Figure  5.1),  or  that  they 
consist  of  straight  line  segments  whose  orientations  with  respect  to  the  ar-axis  are  multiples 
of  7t/4.  We  say  a  routing  is  valid  (for  a  particular  set  of  design  rules)  if  it  satisfies  the  design 
rules. 

We  are  concerned  with  the  following  three  placement  problems: 

1.  Minimum  Separation  Problem.  Given  a  set  of  design  rules  and  a  fixed  offset,  find  the 
minimum  separation  that  allows  a  valid  routing. 

2.  Offset  Range  Problem.  Given  a  set  of  design  rules  and  a  fixed  separation,  find  all  offsets 
that  allow  a  valid  routing. 

3.  Optimal  Offset  Problem.  Given  a  set  of  design  rules,  find  an  offset  minimizing  the 
separation  necessary  for  a  valid  routing. 
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Figure  5.1:  A  rectilinear  routing  connecting  the  modules  P  and  Q. 


Figure  5.2:  The  first  five  barriers  about  pi  when  (a)  arbitrary  wires  are  allowed  and  (b)  only 
rectilinear  wires  are  permitted. 

Before  we  can  relate  these  problems  to  array-searching,  however,  we  need  a  few  more  defini¬ 
tions. 

We  say  that  the  terminals  of  P  and  Q  are  monotone  if  p,  <  q,  for  all  i  or  if  p,  >  q, 
for  all  t.  For  any  P  and  Q,  we  can  assume  that  the  terminals  are  monotone,  since  we 
can  always  partition  the  terminals  into  maximal  monotone  blocks  and  find  a  routing  for 
each  of  the  blocks  independently.  Without  loss  of  generality,  we  assume  p,  <  qt  for  all  i. 
Since  the  terminals  are  monotone,  we  can  assume  that  the  wires  are  also  monotone,  i.e.,  the 
y-coordinate  of  a  wire  is  nonincreasing  as  the  x-coordinate  increases. 

For  1  <  i  <  n  and  1  <  k  <  n  —  i,  we  define  the  k-th  barrier  about  pi  to  be  the  set  of  points 
that  delimit  the  closest  possible  approach  to  pi  of  the  monotone  wire  going  from  p,+*  to  qi+k- 
The  barriers  are  determined  by  the  design  rules.  For  example,  if  the  only  design  constraint 
is  the  lower  bound  on  the  distance  between  wires  (the  case  consider  by  Tompa  in  [Tom80]), 
then  the  barriers  are  composed  of  circular  arcs  and  line  segments,  as  in  Figure  5.2(a).  If,  on 
the  other  hand,  the  design  rules  allow  only  rectilinear  wires,  then  the  barriers  are  rectilinear, 
as  in  Figure  5.2(b). 

Now  consider  the  following  n  x  n  array  A  =  {a.j}.  If  i  <  j,  then  a,j  is  the  height  at  q:  of 
the  ( j  —  z)-th  barrier  about  pi}  as  in  Figure  5.3.  If  i  >  j,  then  at  J  =  0.  For  i  <  j ,  the  value 
aitj  is  the  minimum  separation  possible  for  P  and  Q  given  the  interaction  of  the  terminal  pi 
with  the  wire  running  from  p}  to  qr 
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Figure  5.3:  The  interaction  of  p,  with  the  wire  running  from  pj  to  q3  requires  P  and  Q  be  separated 
by  at  least  dij. 


In  [SD81],  Siegel  and  Dolev  showed  that  the  minimum  separation  possible  for  P  and 
Q  is  simply  the  maximum  entry  in  A.  They  also  showed  that  under  some  very  general 
assumptions  about  barriers,  corresponding  to  a  wide  variety  of  design  rules,  the  array  A 
is  totally  monotone.  (In  fact,  their  proof  is  easily  extended  to  show  A  satisfies  the  Monge 
condition.)  Under  the  assumption  that  each  entry  of  A  can  be  computed  in  constant  time 
(which  is  true  in  most  practical  applications),  Siegel  and  Dolev  then  gave  an  0(n  lg  n)  time 
algorithm  for  the  minimum  separation  problem.  Aggarwal  et  al.  [AKM*87]  reduced  this  to 
0(n)  time  by  simplying  applying  their  linear  time  algorithm  for  computing  the  row  maxima 
of  a  totally  monotone  array. 

In  this  paper,  we  apply  array-searching  to  the  offset  range  and  optimal  offset  problems. 
We  first  show  that  the  offset  range  problem  can  be  solved  in  linear  time  for  most  design 
rules  that  occur  in  practice.  In  particular,  we  consider  design  rules  inducing  what  Siegel  and 
Dolev  refer  to  in  [SD81]  as  a  family  of  similar  concave  barriers.  Such  a  family  of  barriers  is 
characterized  by  a  function  ft(x),  where  for  0  <  x  <  1,  0  <  h(x)  <  1  and  h(x )  is  concave, 
and  for  x  >  1,  h(x)  =  0.  The  height  of  the  k- th  barrier  about  a  terminal  p,  at  a  horizontal 
distance  of  x  from  p,  is  then  given  by  kh(x/k). 

Given  a  fixed  separation  w,  define  the  left  offset  function  L(i,j )  and  the  right  offset 
function  R{i,j)  as  follows:  for  1  <  i  <  n  and  1  <  j  <  n, 


and 


L(i, 


R(i, 


Pi  -  9j  +  {j  -  i)h  MWO'-O) 

if  j  —  t  >  w 

—  oo 

if  j  —  i  <w 

Pi  -  q,  -  {i  -  ;>-1(W(*  -  j )) 

if  f  —  j  >  w 

oo 

if  i  —  j  <  w 

Siegel  and  Dolev  showed  that  the  offset  range  for  Q  is 

nax  L{i,j),  min  R(i,j)\ 

:t,j<n  1<i,  ]<n  J 


They  also  proved  the  following  theorem. 
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Theorem  5.1  ([SD81])  Suppose  r  >  0,  s  >  0,  and  j  >  i  +  r.  If  L(i,j)  <  L(i  +  r,  j),  then 
L(i,j  +  s)  <  L{i  +  r,j  +  s).  Similarly,  ifR{i,j)  >  R(i  +  r,j),  then  R(i,j+$)  >  R(i+r,j  +  s). 

Now  let  L  =  {lij}  denote  the  n  x  n  array  where  l,tJ  =  L(j,i)  and  let  R  =  {r,-j}  denote 
the  nxn  array  where  =  R(j,  i).  An  immediate  corollary  of  Theorem  5.1  is  the  following. 

Corollary  5.2  The  arrays  L  and  R  are  totally  monotone.  Furthermore,  if  L(i,j)  and  R(i,j ) 
can  be  computed  in  constant  time,  then  the  offset  range  problem  can  be  solved  in  linear  time 
by  computing  the  row  maxima  of  L  and  the  row  minima  of  R. 

Proof  Consider  any  2x2  minor  of  L ,  corresponding  to  rows  i  and  k,  i  <  k,  and  columns 
j  and  /,  j  <  l.  If  l  >  i,  then  U,i  =  L(l,i)  =  —  oo,  which  implies  the  minor  is  monotone.  If,  on 
the  other  hand,  /  <  i,  then  by  Theorem  5.1,  Z,j  <  /,,/  implies  lkj  <  lkj,  which  again  means 
the  minor  is  monotone.  Thus,  L  is  totally  monotone.  We  can  show  R  is  totally  monotone 
in  a  similar  manner. 

If  we  can  compute  any  entry  of  L  or  R  in  constant  time,  then  we  can  apply  the  results  of 
[AKM*87]  and  compute  L' s  row  maxima  and  R' s  row  minima  in  linear  time.  We  can  then 
compute  the  maximum  of  L's  row  maxima  and  the  minimum  of  R's  row  minima  in  linear 
time  to  obtain  the  offset  range  for  Q.  ■ 

This  corollary  generalizes  the  results  of  Siegel  [SieSS],  who  gave  an  0(n )  time  algorithm 
for  the  offset  range  problem  when  barriers  are  polygonal  in  nature. 

Passing  now  to  the  optimal  offset  problem,  Siegel  and  Dolev  showed  in  [SD81]  that  this 
problem  can  be  solved  for  a  particular  set  of  design  rules  in  lg  n )  time  if  the  minimum 

separation  problem  can  be  solved  for  these  design  rules  in  0(^(n))  time.  Consequently, 
Corollary  5.2  gives  us  an  0(n  lg  n)  time  algorithm  for  the  optimal  offset  problem  for  most 
design  rules  used  in  practice.  Though  it  remains  open  whether  an  o(nlgn)  time  algorithm 
for  this  problem  can  be  obtained  for  the  class  of  design  rules  we  consider,  Mirazaian  [Mir87] 
has  provided  an  elegant  0(n)  time  algorithm  for  the  case  of  rectilinear  wiring;  however,  his 
algorithm  exploits  the  properties  associated  with  such  a  wiring,  and  it  is  unlikely  that  his 
techniques  can  be  extended  to  other  wiring  models. 

6  Open  Problems 

In  this  paper,  we  present  a  framework  that  allows  us  to  obtain  efficient  algorithms  for  wide 
variety  of  problems  involving  quadrangle  inequalities.  We  leave  a  large  number  of  problems 
unresolved,  however.  These  open  problems  may  be  divided  into  two  types:  those  relating 
directly  to  array-searching  and  those  involving  applications  of  array-searching. 

Among  the  open  problem  relating  directly  to  array-searching  are  the  following: 

1.  In  Subsection  2.1,  we  presented  an  O(nlgn)  time  algorithm  for  computing  the  plane 
maxima  of  an  n  x  n  x  n  three-dimensional  totally  monotone  array.  This  is  the  best 
upper  bound  we  have  on  the  time  necessary  to  solve  this  problem,  even  for  the  special 
case  of  Monge-composite  arrays,  whereas  the  only  lower  bound  known  for  either  of 
these  problems  is  fi(n).  Note  that  our  algorithms  for  computing  the  plane  maxima 
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of  a  three-dimensional  array  A  only  make  use  of  the  monotonicity  of  A  and  the  total 
monotonicity  of  its  planes.  It  is  conceivable  that  an  algorithm  that  takes  advantage  of 
the  stronger  Monge  condition  might  be  able  to  solve  the  plane  maxima  problem  for  A 
in  o(nlgn)  time. 

2.  In  a  similar  vein,  the  only  upper  bound  known  for  the  time  necessary  to  compute  the 
plane  maxima  of  an  n  x  n  x  •  •  •  x  n  d-dimensional  totally  monotone  array  is  0(n  \gd~ 2  n), 
whereas  Q(dn)  remains  the  only  lower  bound  known  for  this  problem.  Reducing  this 
upper  bound  to  something  polynomial  in  d  would  be  a  major  improvement.  Again,  it 
is  conceivable  that  the  stronger  Monge  condition  might  be  of  use  in  this  regard. 

3.  In  [AK88],  Aggarwal  and  Klawe  defined  an  n  x  n  staircase-monotone  array  to  be  a 
totally  monotone  array  with  the  entries  in  some  contiguous  region  of  the  array  replaced 
by  —  oo’s,  so  that  the  array  is  no  longer  totally  monotone.  (This  region  is  constrained 
to  lie  in  the  upper  left  corner  of  the  array  and  to  be  bounded  by  a  staircase  growing 
up  and  to  the  right.)  Aggarwal  and  Klawe  showed  that  the  row  maxima  problem  for 
such  an  array  can  be  solved  in  O(nlglgn)  time  and  gave  a  number  of  applications  in 
which  such  arrays  arise.  Klawe  and  Kleitman  [KK88]  improved  the  time  bound  on  this 
problem  to  0(nct(n)),  where  a(n)  is  the  inverse  Ackermann  function.  The  only  lower 
bound  known  for  this  problem  is  fi(n). 

As  for  the  applications  of  array-searching  we  present  in  Sections  3-5,  we  list  the  open 
problems  in  these  areas  in  their  respective  sections.  Clearly,  reducing  any  of  the  upper  bounds 
given  for  the  open  problems  listed  above  will  result  in  the  improvement  of  the  corresponding 
application  algorithms,  but  it  is  quite  possible  that  exploiting  additional  properties  asso¬ 
ciated  with  the  individual  application  problems  may  help  in  reducing  their  computational 
complexity  as  well. 

In  addition  to  the  applications  that  we  describe  in  this  paper,  there  are  also  a  number 
of  problems  that  seem  closely  related  to  the  notion  of  array-searching.  Recall  that  Frances 
Yao’s  dynamic  programming  algorithm  requires  that  the  weight  function  w{i,j)  satisfy  the 
quadrangle  inequality  and  an  additional  monotonicity  constraint  (namely,  w(i',j')  <  w(i,j ) 
when  i  <  i‘  <  j‘  <  j).  Our  algorithm  for  Yao’s  problem,  given  in  Subsection  4.1,  also 
requires  that  the  weight  function  satisfy  this  monotonicity  constraint.  This  constraint  seems 
somewhat  restrictive;  there  are  a  number  of  problems,  quite  similar  to  Yao’s  dynamic  pro¬ 
gramming  problem,  that  involve  weight  functions  satisfying  the  quadrangle  inequality  but 
not  the  additional  monotonicity  constraint.  For  example,  Gilbert  [Gil79]  has  shown  that  a 
minimum  weight  triangulation  of  a  convex  n-gon  can  be  found  in  0(n3)  time  using  a  dy¬ 
namic  programming  algorithm.  Though  the  lengths  of  the  diagonals  added  in  triangulating 
a  convex  n-gon  do  satisfy  the  quadrangle  inequality,  they  do  not  obey  Yao’s  monotonicity 
constraint,  thus  neither  Yao’s  techniques  nor  our  own  are  directly  applicable.  Furthermore, 
by  incorporating  the  notion  of  optimal  t- ary  trees,  Gilbert’s  algorithm  can  be  extended  to  ob¬ 
tain  an  0(n 3  lg  m)  time  algorithm  for  the  problem  of  partitioning  a  convex  n-gon  into  convex 
m-gons  such  that  the  sum  of  the  perimeters  of  these  m-gons  is  minimized  [LLS87].  However, 
again,  despite  the  similarity  with  Yao’s  dynamic  programming  problem,  no  o(n3Igm)  time 
algorithm  for  this  problem  is  known. 
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